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ABSTRACT 

We present results of a systematic study of failing core-collapse supernovae and the formation of stellar-mass 
black holes (BHs). Using our open-source general-relativistic 1.5D code GRID equipped with a three-species 
neutrino leakage/heating scheme and over 100 presupernova models, we study the effects of the choice of 
nuclear equation of state (EOS), zero-age main sequence (ZAMS) mass and metallicity, rotation, and mass- 
loss prescription on BH formation. We find that the outcome, for a given EOS, can be estimated, to first order, 
by a single parameter, the compactness of the stellar core at bounce. By comparing protoneutron star (PNS) 
structure at the onset of gravitational instability with solutions of the Tolman-Oppenheimer-Volkof equations, 
we find that thermal pressure support in the outer PNS core is responsible for raising the maximum PNS 
mass by up to 25% above the cold NS value. By artificially increasing neutrino heating, we find the critical 
neutrino heating efficiency required for exploding a given progenitor structure and connect these findings with 
ZAMS conditions, establishing, albeit approximately, for the first time based on actual collapse simulations, the 
mapping between ZAMS parameters and the outcome of core collapse. We also study the effect of progenitor 
rotation and find that the dimensionless spin of nascent BHs may be robustly limited below a* = Jc/GM 2 = 1 
by the appearance of nonaxisymmetric rotational instabilities. 



Subject headings: black hole physics - equation of state - hydrodynamics 
mass-loss - stars: neutron - stars: supernovae: general 



neutrinos - stars: evolution - stars: 



1. INTRODUCTION 

Massive stars with zero-age main sequence (ZAMS) masses 
M Z ams in the range of 8- 10M Q < M ZA ms < 100- 150M Q 
end their lives with the gravitationally induced catastrophic 
collapse of their electron-degenerate iron core to nuclear 
densities. There, the nuclear equation of state (EOS) stiff- 
ens and stabilizes the inner core, which overshoots its new 
equilibrium and bounces back, launching a hydrodynamic 
shock. The shock initially races through the still collaps- 
ing outer core, but soon stalls and turns into an accretion 
shock (at r ~ 100 -200 km) due to the dissociation of heavy 
nuclei at the s hock front an d neutrino losses from the post- 
shock region ([BetheJ 119901) . The shock must be revived 
to drive a core-collapse supernova (CCSN) and the pre- 
cise nature of the responsible CCSN mechan ism has been 
a topic of intense research for decades ( e.g.. lArnetg 119661: 
Colgate & White! [19661 bethe & Wilson! 119851: Uanka et al l 
2007tlBurrows et al.ll2006l l2007bt iMurphv & Burrowsll2008l: 
Marek & Jankal 120091: iNordhaus et al.1 120101 and references 
therein). 

A neutron star (NS) is left behind by a CCSN that explodes 
soon after bounce and successfully unbinds its stellar man- 
tle. However, a stellar-mass black hole (BH) may be the out- 
come: (1) if in a successful, but perhaps weak, CCSN fall- 
back accretion pushes the nascent NS over its mass limit, (2) 
if nuclear phase transitions during protoneutron star (PNS) 
cooling occur or if PNS cooling reduces pressure support in 
a hyper-massive PNS, or (3) if the CCSN mechanism lacks 
efficacy and fails to revive the shock and continued accretion 
pushes the PNS over its maximum mass. In this last chan- 
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nel to a stellar-mass BH, there is no electromagnetic (EM) 
signal other than the disappearance of the original star. It is 
such "unnovae" dKochanek et al.ll2008l) . failing CCSNe, that 
are the topic of this paper. 

In ordinary massive stars that hydrostatically form degener- 
ate ir on cores, BH formation, in any scenari o, is never prompt 
(e.g. jBurrowsl[T^lolt^O , Connodl20T0l) . It is always pre- 
ceded by an extended P NS phase giving rise to copious emis- 
sion of both neutrinos dBurrowsl M988T iBeacom et al.ll200ft) 
and gravitational waves (lOttl 120091) until the PNS is engulfed 
by the BH horizon. The EM silence expected in a failed 
CCSN may be broken after all, if sufficient and appropriately 
distributed angular momentum is present to allow for a Keple- 
rian accretion d isk to form near the BH, permitting a collapsar 
(iWooslevI 19931) gamma-ray burst (GRB) central engine to op- 
erate and drive relativistic outflows. 

It is currently unclear what fraction of massive stars 
form BHs and through which channel. Preexplosion ob- 
servations of progenitors of su ccessful CCSNe su ggest pro- 
genitor masses < 17-20M Q (ISmartt et all 120091) for stan- 
dard Type I I -P sup ernovae. Assuming, as suggested by 
ISmartt et all (120091) . that most other CCSNe fail or make 
BHs after a successful explosion, this would correspond 
to a BH fraction of < 30 %-35 % of massive stars above 
8Mq. However, al ternative interpret ations exist and have 
been summarized bylSmith et al.l (120101). Theoret i cal wo rk by 
i Timmes et all (119961). iFrv erJ (1 19991). i Heger et all (120031) . and 
lEldridge & Tout! (120041) provided rough estimates on the out- 
comes of stellar collapse as a function of progenitor ZAMS 
mass and metallicity. Leavin g effects due to binary evolu- 
tion aside, IZhang et all (120081) performed an extensive study 
of fallback in artificially driven spherically symmetric CCSN 
explosions and estimated that zero-metallicity stars form BHs 
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in 20%-50% of all core-collapse events with an average 
BH mass of 6- 10M Q . For solar-met allicity stars, due to 
increased mass loss during evolution, IZhang et~aT] d2008l) 
found BHs to form at a significantly lower rate and initial 
mass. They predict BH fractions in the range of 10%-25% 
with typical initial BH masses of 3M . This is in rough 
agreement with previous pop ulation synthesis c alcula tions of 
iFrver & Kalogeral (12001b and lBelczvnski et al.1 (120021) . 

Early spherically symmetric (one-dimensional, ID) simu- 
lat ions of BH for mati on in failing CCSNe were c arried out 
by iWilsonl (11971b and Ivan Riper & Arnettl dl978l) . iBurrowsl 
119881) performed a set of quasi-stationary ID PNS accre- 
tion and cooling simulations to investigate the possibility of 
BH formation in SN 1987 A. Delayed BH formation (by tens 
of seconds), du e to, e.g., a nuc l ear EOS phase transition, 
was studied by Baumga rte et al.l dl996albl) . More recently, 
ID full Boltzmann neutrino radiation-hydrodynamics calcu- 
lations of failing nonrot ating CCSNe wer e carried out by 
iLiebendorfer et alj (120041) and more recentlv | Sumivoshi et alJ 
d2007L I2008l " l2009h and iFischer et alJ d2009j). These stud- 
ies provided detailed neutrino signature predictions for BH- 
forming core-collapse events. However, owing to the 
complexity and computational expense of such Boltzmann- 
transport calculations, these groups could consider only very 
limited sets of progenitor models and EOS. Simplified ax- 
isymmetric (2D) simulations of BH formation in rotating 
core collapse were first performe d in a series of papers by 
ISekiguchi & Shibatal d200H l200l) and IShibata & Sekiguchil 
d2005l) . These authors used simplified EOS, no neutrino 
treatment, and artificially constructed initial conditions and 
found prompt BH formation. Recently, the same authors 
performed a small set of 2D simulations with a finite- 
temperature nuclear EOS an d a leakage scheme for neutrinos 
dSekiguchi & Shibatal2010l) and considered collapse, BH for- 
mation, and subsequent evolution in an artificially constructed 
progenitor with an iron core mass of ~ 13M Q and constant 
specific entropy of 8 ^/baryon, initial conditions that are in- 
consistent with those at the precollapse stage of CCSN pro- 
genitors. 

In this paper, our focus is on studying and establishing the 
systematics of failing CCSNe and BH formation. For this, we 
employ the spherically symmetric general-rela tivistic (GR) 
open-source code GRID dO' Connor & Ottl2010l) that can han- 
dle rotation in an approximate angle-averaged way ("1.5D") 
and sacrifice accuracy in the neutrino treatment by employ- 
ing an efficient energy-averaged three-species neutrino leak- 
age/heating scheme instead of full transport. The efficiency of 
GRID enables us to perform more than ~ 700 collapse calcu- 
lations, investigating for the first time in detail the effects of 
variations in nuclear EOS, progenitor ZAMS mass and metal- 
licity, neutrino heating efficiency, and precollapse rotational 
configuration. We employ four different finite-temperature 
nuclear EOS and draw a total of 106 progenitor models from 
seven stellar evolution studies. 

In Section|2] we review the features of our 1.5D GR hydro- 
dynamics code GRID, discuss our neutrino leakage/heating 
scheme, and introduce the set of employed EOS. Section[3]in- 
troduces our progenitor model set, nu meric al grid setup, and 
precollapse rotational setup. In Section l4~Tl we introduce key 
aspects of failing CCSNe and BH formation by discussing the 
evolution of BH formation in a fiducial nonrot ating 40M Q 
solar-metallicity progenitor. We go on in Section l-Ol to study 



the influence of the EOS and thermal effects on the time to 
BH formation and on the maximum (baryonic and gravita- 
tional) mass of the PNS. We discover that for nuclear EOS 
with physically plausible stiffness, the maximum (baryonic 
and gravitational) mass of the PNS is always greater than the 
corresponding cold NS mass and discuss that the difference is 
due entirely to thermal pressure support of material in the hot 
outer PNS core. This effect is strongest for the softest con- 
sidered EOS and decreases with increasing EOS stiffness. In 
Section 14.41 we analyze the impact of variations in progeni- 
tor structure on the time to BH formation and the maximum 
PNS mass in failing CCSNe. We find that the postbounce dy- 
namics can be predicted rather robustly by a single parameter, 
the compactness of the progenitor structure at core bounce. 
The same appr oximate single-parameter dependence emerges 
in Section 1431 where we determine the neutrino heating effi- 
ciencies required (modulo ignored multi-dimensional effects) 
to induce a neutrino-driven explosion in a larg e s et of pro- 
genitors. The combined results of Sections [4.41 and 1431 allow 
us to make predictions on the outcome of core collapse for 
pro genit ors with varying ZAMS mass and metallicity in Sec- 
tion [4j6] As we discuss in that section, mass loss may be the 
greatest uncertainty in conn ectin g ZAMS parameters to core- 
collapse results. In Section |4~71 we present results from the 
first rotating BH formation simulations in the CCSN context. 
Varying the precollapse rotation rate in Section l4.7.1l we find 
that, not unexpectedly, increased rotation leads to a delay of 
BH formation and greater maximum PNS masses. We also 
observe that the birth spin of Kerr BHs in nature appears to be 
robustly limited to values below a* = J/M 2 < 0.9 by the likely 
appearance of nonaxisymmetric rotational instabilities that re- 
distribute or radiate angular momentum. This finding re quires 
confirmation by 3D simulations. We go on in Section 14.7.21 
to discuss the collapse evolution of a set of progenitors that 
were evolved from the ZAMS with a 1.5D treatment of ro- 
tation and discuss their viability as collapsar-type long-GRB 
progenitors. Finally, in Section|5] we critically summarize our 
work and conclude. 

2. METHODS 

2.1. GR Hydrodynamics 

GRID dO' Connor & Ottll2010l) is a spherically symmetric 
GR hydrodynamics code developed for the study of stellar 
collapse and BH formation. It is available for download 
at |http : //stellarcollapse . org , GR ID, based on 
the pre vious work of IGourgoulhonl d 199 11) and lRomero et alJ 
d!996l) . is Eulerian and uses the radial gauge - polar slicing 
coordinates that have the simplifying property of a vanishing 
shift vector. Here, we briefly outline the basics of the cur- 
vature and hydrodynam ics equations and refer the reader to 
IP' Connor &"Ottl (l2010h for full details and derivation. We 
assume spacelike signature (-,+,+,+) and, unless noted oth- 
erwise, use units of G = c = M & = 1. The metric of GRID is 
given by the line element 

ds 2 = -air, tfdt 2 + X(r, tfdr 2 + r 2 dfl 2 , ( 1 ) 

where a{r, t) = exp($(r, t)) with <£>(r, t) being the metric poten- 
tial and X(r,t) = [1 -2m(r)/r]" 1//2 , where m(r) is the enclosed 
gravitational mass. We assume an ideal fluid with stress en- 
ergy given by 

T» v ^phuW+g^P, (2) 
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where p is the matter density, P is the fluid pressure and 
h= 1+e + f/pis the specific enthalpy with e being the spe- 
cific internal energy. m m in Equation (f2]i is the 4- velocity of the 
fluid, and without rotation, taken to be = (W/a, Wv r ,Q,Q), 
where W = [1 — v 2 ] -1 / 2 is the Lorentz factor and v = Xv r is 
the physical velocity. For a given matter configuration, the 
Hamiltonian and momentum constraint equations give differ- 
ential equations for both m(r) and $(r), 



m(r) ■ 



:4tt / {phW z 
'o 



-P+T»y 2 dr', 



(3) 



m(r',t) 



+ 4irr'(phX 2 u r 'u r ' +P+t£) 



dr'+$ , 
(4) 

where $o is determined by matching the metric at the star's 
surface to the Schwarzschild metric. The neutrino terms, t£ 
and r^, ac count for trapped neutrin os and their detailed form 
is given in |Q' Connor & Ottl (I2010I) . We obtain the fluid evo- 
lution equations by expanding the local fluid rest-frame con- 
servation laws, V^T^ = and V^/^ = 0, in the coordinates 
of GRID. The conservation laws become 



X 



■s 



(5) 



where U = (D, DY e , S r , t) are the conserved variables, given 
in terms of the primitive fluid variables p, Y e , e, P, and v as 




XpW 
XpWY e 
phW'v 
phW 2 -P-D 



(6) 



The spatial fluxes in Equation §5^ are given by, 
F = (Dv, DY e v, S'v+P, S r -Dv), 
and the source terms 



(7) 
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0, Rl, (S'v-T-D)aX(&TrrP+^) + aPX^ 



Q^+Qr' 
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where the Rs and Qs are neutrino sources and sinks which 
arise from the neu trino leakage scheme a nd neutrino pressure 
contributions (see lO'Connor & Ottll2010l for details). 

The evolution equations (Equation |5]l are firs t spatially dis- 
critized using a fi nite-volume scheme (e.g., iRomero et al.l 
[l996t iFontl [20081). The piecewise parabolic method 
(IColella & Woodward! [1981 is used to reconstruct the state 
variables at the cell interfaces and the HLLE Riemann solver 
(lEinfeldtll 19881) is employed to determine the physical fluxes 
through these interfaces. The evolution equ ations are inte - 
grated forward in time via the method of lines (iHvmanl 19761) . 
using standard second order Runge-Kutta time integration 
with a Courant factor of 0.5. After updating the conserved 
variables, they are inverted via a Newton-Raphson scheme to 
obtain the new fluid state variables. 

In spherical symmetry, rotation can be included by assum- 
ing constant angular velocity fi on spherical shells (shellular 
rotation) and including an angle-averaged centrifugal force 



in the radial momentum equation. This is commo n prac- 
tice in stellar evolution codes (e.g.. lHeger et al. 2000) an has 
also b een applied to Newtonian ID stellar collaps e calcula- 
tions dAkEamaetiD Exjl EompsonetiD|2005l). We in- 
clude a GR variant of this "1.5D" rotation treatment in GRID 
(1) by adding an evolution equation for the generalized spe- 
cific angular momentum 5^ = phW 2 v tf> r, (2) by including an 
effective centrifugal force in the equation for S r , and (3) by 
modifying the expressions for the 4-velocity, the Lorentz fac- 
tor, and the differential equation for the metric potential to 
account for rotation. Full details as well as a demonstra- 
tion of conserva t ion of angular momentum can be found in 
lO'Connor & Ottl (12010I). Note t hat, a s may be expected and 
was demonstrated by loit etail d2006h . the 1.5D approxima- 
tion becomes less accurate with increasing spin and quantita- 
tive results are reliable only for low rotation rates. 

2.2. Neutrino Treatment 

Neutrino effects are crucial in stellar collapse and 
should ideally be included via a comput ationally expensive 
GR B oltzmann transport treatment (e.g., ILiebendorfer et all 
120041) . However, since our aim is to perform an extensive 
parameter study with hundreds of simulations, we choose to 
resort to a less accurate, but much more computationally effi- 
cient leakage and approximat e heating scheme for neu trinos. 
Details of this are laid out in lO'Connor & Ottl (120101) . Here 
we review only its most salient features. 

Before core bounce, neutrinos deleptonize the collaps- 
ing core, reducing the electron fraction Y e an d, as a conse - 
quence, the size of t he homologous inner core (lBethelll990l) . 
ILiebendorfer! (120051) showed that the prebounce Y e can be pa- 
rameterized as a function of density and that this parame- 
terization varies little between progenitor stars. We follow 
this prescription for prebounce deleptonization and use the 
Y e (p) fit parameters of his G15 model. Following bounce, 
this simple parameterization becomes inaccurate and cannot 
capture the effects of neutrino cooling, deleptonization, and 
neutrino heating. Hence, we switch to a leakage scheme that 
uses elements of what was la id out by iRuffert et al.l (119961) 
and lRosswog & Li ebendorferl (120031) . We consider three neu- 
trino species, v e , D e , and v x = {1/^,9^,^,9^. Neutrino pairs 
of all species are made in thermal processes of which we in- 
clude electron-posit ron pair annihilation and plasmon decay 
(IRuffert et al.lll996l) . In addition, charged-current processes 
lead to the emission of v e s and D e s. The leakage scheme pro- 
vides approximate energy and number emission rates that are 
inserted into GRID'S evolution equations via source terms in 
Equation ©, R% , Q^, and Q" T ' E (lO'Connor & Ottll20Toh . 

We include neutrino heating vi a a paramete rized charged- 
current heating scheme based on iJankal (1200 ll) . The heating 
rate at radius r is 



Mr) P Y I 1 

47TT Z TO,, \ t v . 



-2t„. 



(9) 



where /h eat is a scale factor that allows for artificially in- 
creased heating, L Vj {r) is the neutrino luminosity interior to 
r, r Vi is the optical depth, determined through the leakage 
scheme, Cheats is the energy-averaged absorption cross sec- 
tion, and Xj is corresponding mass fraction of the neutrino 
reaction (X p for v e capture on protons and X„ for v e capture 
on neutrons). (1 /F Vj ) is the mean inverse flux factor which we 
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approximate analytically as a function of the optical depth r 
by compari ng to the angle-d ependent radiation transport cal- 
culations of lOtt et all (12008b . 

We include in our simulations the stabilizing effect of 
neutrino pressure in the opticall y thick PNS core via an 
ideal Fermi-gas appro ximation dLiebendorfer et all 12005b 
|Q' Connor & Ottl 120101) . Leaving out this pressure contri- 
bution leads to ~ 5% smaller maximum gravitational PNS 
masses. We also include terms due to neutrino pressure and 
radiation-field energy in the calculation of the gravitational 
mass (Equation [3j and of the metric potential (Equation HJ. 
Since our leakage scheme does not treat neutrino energy sep- 
arately from the internal energy of the fluid, including the en- 
ergy of the neutrino gas in the former equations is not fully 
consistent with our present approach. This error was dis- 
covered and corrected after all simulations were performed. 
However, a set of test calculations showed that the error leads 
to an underestimate of the maximum gravitational PNS mass 
of only ~ 2% which is well wit hin the error of the overall 
leakage scheme (see also Section l4~2l . 

2.3. Equations of State and Maximum Neutron Star Masses 

We include multiple finite-temperature nuclear EOS in this 
study to explore the dependence of postbounce evolution and 
BH fo rmation on EOS propertie s. The Lattimer-Swesty (LS) 
EOS dLattimer & Swestvl 1 199 lb is based on the compress- 
ible liquid-droplet model, assumes a nuclear symmetry en- 
ergy S v of 29.3 Me V, and comes in three variants with differ- 
ent values of the nuclear incompressibility of K s = 180MeV 
(L S 1 80), 220 Me V (L S220), and 375 MeV (LS375). The EOS 
of IShen et al.l (Il998allbh (HShen EOS), on the other hand, is 
based on a relativistic mean-field model, has S v = 36.9 MeV 
and K s = 281 MeV. More details on these EOS and t h eir im - 
plementation in GRID is given in |Q' Connor & Ottl d2010l) . 
The EOS tables and driver routines employed in this study 
are available for download at http://stellarcollapse.org 

By solving the Tolman-Oppenheim er-Volkoff (TOV) equa- 
tions (lOppenheimer & Volkoffl[l939l) with T = 0.1 MeV and 
assuming neutrinoless [3 equilibrium we determine the neu- 
tron star baryonic and gravitational mass-radius relationships 
that each of these four EOS produces and that are depicted by 
FigureQ] The maximum gravitational (baryonic) neutron star 
masses are - 1.83 M (~ 2.13M ), ~ 2.O4M (~2.41M Q ), 

- 2.72M Q (~ 3.35M Q ), and ~ 2.24 M (~ 2.61 M ) for 
LSI 80, LS220, LS375, and HShen, respectively. The co- 
ordinate radii of these maximum-mass stars are — 10.1km, 
~ 10.6km, — 12.3km and — 12.6km, respectively. 

The above maximum neutron star masses hold only for non- 
rotating cold NSs. As we will discuss in detail in Section 
14.31 the PNSs at the heart of the failing CCSNe considered in 
this work are much hotter. They have central temperatures of 

— 1 0- 20 MeV and tens of MeV in their outer core and mantle. 
Thermal effects have a significant effect on their maximum 
masses. 

In this study, we do not consider hypero nic EOS, e.g.. 
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the hy peronic extension of the HShen EOS by llshizuka et al.l 
(120081) . or EOS in volving other phases o f nuclear matter, e.g., 
quarks and pions [Nakazato et al.l ( 120101) . Such EOS are po- 
tentially interesting in failing CCSNe, since their exotic com- 
ponents lead to a softening of the EOS at high density, p oten- 
tially accelerating BH formation (ISumivoshi et al.ll2009l) . We 
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Figure 1. Baryonic (left) and gravitational (right) neutron mass-radius re- 
lations for various hot nuclear EOS. The temperature is taken to be constant 
throughout the star at T = 0. 1 MeV and the electron fraction is determined 
through neutrinoless /3-equilibrium with an imposed minimum of 0.05 due to 
table constraints. 

also do not consider EOS that include QCD phase transitions 
that too may lead to early PNS co llapse and potenti ally to a 
second bounce and neutrino burst (Sagert et al. 2009). 

3. MODEL SETUP 

3.1. Pre supernova Data 

We make use of single-star nonrotating presupernova mod- 
els fro m several stella r evolution stud ies : IWooslev & Weaverl 
(fl995l) (WW 95) iWooslev et al.l (120021) (WHW02) 
iLimongi & Chief! (120061) (LC06A/B) and IWooslev & Hegeri 
(120071) (WH07). Each of these studies evolved stars with a 
range of ZAMS masses at solar metallicity (Z , hereafter 
denoted with prefix s in model names) up until the onset 
of core collapse. In addition to solar metallicity, WHW02 
evolved stars with ultra low metallicity, 1O~ 4 Z (denoted 
by prefix u) and zero metallicity (denoted by prefix z). 
Rotation is of relevance in stellar evolution and stellar 
evolutionary processes affect the rotational configuration at 
the presupernova stage. In order to study BH formation, BH 
birth properties and their impact on a potential subsequent 
evolution to a GRB in such spinning progenitor s, we draw 
representa tive models from Heger et ail (120001) (HLW00) 
and from IWooslev & Hegeri d2006j) (WH06) who included 
rotation in essentially the same way as we do in GRID. 

In Table [T] we list key parameters for all models in our set. 
These include presupernova mass, iron core mass (which we 
define as the baryonic mass interior to Y e = 0.495), and the 
bounce compactness £2.5. The latter is defined as 



£m - 



M/M Cl 



R(M hm =M)/ 1000km 



(10) 



where we set M = 2.5M . R{M hmy = 2.5M ) is the radial 
coordinate that encloses 2.5-M at the time of core bounce. 
£2.5 gives a measure of a progenitor's compactness at bounce. 
We choose M = 2.5M as this is the relevant mass sc ale f or 
BH formation. £2.5 is, as we shall discuss in Section l4~4l a 
dimensionless variable that allows robust predictions on the 
postbounce dynamics and the evolution of the model toward 



BH formation. The evaluation of £2.5 at core bounce is cru- 
cial, since this is the only physical and unambiguous point in 
core collapse at which one can define a zero of time and can 
describe the true initial conditions for postbounce evolution. 
Computing the same quantity at the precollapse stage leads to 
ambiguous results, since progenitors come out of stellar evo- 
lution codes in more or less collapsed states. Collapse washes 
out these initial conditions and removes ambiguities. 

We point out (as is obvious from Table [U that there 
is a clear correlation between iron core mass and bounce 
compactness. Since the effectiv e Chandrasekhar mass in- 
creases due to thermal corre ctions dBurrows & Lattimerl 19831: 
iBaron & Coopersteinll990t) . more massive iron cores are hot- 
ter. Hence, progenitors with greater bounce compactness re- 
sult in higher- temperature PNSs. 

One of the most uncertain, yet most important, vari- 
ables in the evolution of massive stars is the mass-loss 
rate. Mass loss can vary significantly over the life of 
a star. Current estimates of mass loss, either theoreti- 
cal or based on fits to observational data, can depend on 
many parameters, including mass, radius, stellar luminos 
ity, effective surface temperature, surfa ce hydrogen and he 



Table 1 

Initial Models 



lium abundance, and stellar me tallicitv (Ide Jager et al 



Nieuwenhuiizen & de JageJ[T9 90; Well stein & La nger 



1988 



1999 



Nugis & Lamersl l2000t IVink & de Koteri 120051) . The mass- 
loss rate is uncertain in both the massive O-star and in the 
stripped-envelope Wolf-Rayet (W-R) star stage. O-star winds 
are expected to be responsible for the partial or complete re- 
moval of the hydrogen envelopes of massive stars. Recent ob- 
servational results suggest that the rates used in current stel- 
lar evolution models may be too high by factors of 3 — 10 if 
clumped winds are considered correctl y dBouret et al.l 120051 
iFullerton et all l2006t iPuls et alJ l2006h . With the reduced 
rates, WR stars would be difficult to make in standard single- 
star evolut ion and would req uire binary or eruptive mass-loss 
scenarios dSmithetal.1120101) . 

In Figure [2] we plot the mass-loss-induced mapping be- 
tween ZAMS mass and presupernova mass for the ensemble 
of nonrotating progenitors listed in TableQ] WW95 models do 
not include mass loss - the presupernova models of this study 
have a mass equal to the ZAMS mass. WHW02 and WH07 
emplo y the mass-loss rates of Nieuwenhuijzen & de Jagerl 
(119901) and IWellstein & Langerf d 19991) and use significantly 
reduced rates for low and zero metallicity stars. The u 
and z models of WHW02 have almost no mass loss and 
their presupernova masses are very close to their ZAMS val- 
ues. The solar-metallicity stars of the sWHW02 and sWH07 
model sets have significant mass loss, generally scaling with 
ZAMS mass. The most massive stars in these model sets 
have presupernova masses that are a small fraction of the 
initial ZAMS mas s. For main sequence and giant phases, 
Limongi & Chieffil (12006ft a dopt mass-los s rates following 



Vink et al. I (I2000L 120011) and Ide Jager eT al. (1988) . For W-R 
stars, they either use the mass-loss rates of iNugis & Lamers 



( 200C ) (hereinafter referred to as LC06 A models) or lLangei 
(119891) (LC06B models). The latter are close to the values 
used for solar-metallicity stars in the WHW02 and WH07 
model sets. The difference in the LC06A and LC06B mass- 
loss rates is roughly a factor of two. This, as portrayed by 
Figure [2] and evident from Table Q] can significantly alter the 
total mass at the onset of collapse and also has a strong effect 
on the iron core mass and bounce compactness. 



1V1UUC1 


M ZAMS 


M pre-SN 


lv± te core 


6.5" 




(W ) 


(M Q ) 




s20WW95 


20 


20.0 


1.74 


0.383 


S25WW95 


25 


25.0 


1.77 


0.416 


S40WW95 


40 


40.0 


1.98 


0.583 


S15WHWU2 


15 


12.6 


1.55 


0.150 


S20WHW02 


20 


14.7 


1.46 


0.127 


S25WHW02 


25 


12.5 


1.62 


0.326 


S30WHW02 


30 


12.2 


1.46 


0.223 


s35WHW02 


35 


10.6 


1.49 


0.205 


S40WHW02 


40 


8.75 


1.56 


0.266 


S75WHW02 


75 


6.36 


1.48 


0.112 


U20WHW02 


20 


20.0 


1.57 


0.338 


U25WHW02 


25 


25.0 


1.53 


0.223 


u30WHW02 


30 


30.0 


1.58 


0.326 


U35WHW02 


35 


35.0 


1.85 


0.664 


u40WHW02 


40 


40.0 


1.90 


0.719 


u45WHW02 


45 


44.9 


1.96 


0.655 


U50WHW02 


50 


49.8 


1.83 


0.574 


u60WHW02 


60 


59.6 


1.88 


0.623 


U75WHW02 


75 


74.1 


2.03 


1.146 


Z20WHW02 


20 


20.0 


1.48 


0.163 


Z25WHW02 


25 


25.0 


1.81 


0.404 


Z30WHW02 


30 


30.0 


1.50 


0.221 


Z35WHW02 


35 


35.0 


1.79 


0.560 


Z40WHW02 


40 


40.0 


1.90 


0.720 


S25LC06A 


25 


16.2 


1.43 


0.204 


S30LC06A 


30 


12.8 


1.48 


0.274 


S35LC06A 


35 


11.8 


1.48 


0.242 


S40LC06A 


40 


12.4 


1.50 


0.339 


s60LC06A 


60 


16.9 


1.63 


0.603 


S80LC06A 


SO 


22.4 


1.67 


0.628 


S120LC06A 


120 


30.5 


1.91 


0.905 


S40LC06B 


40 


6.82 


1.51 


0.322 


s60LC06B 


60 


5.95 


1.35 


0.163 


S80LC06B 


80 


6.04 


1.46 


0.185 


S120LC06B 


120 


6.12 


1.24 


0.143 


S20WH07 


20 


15.8 


1.55 


0.288 


S25WH07 


25 


15.8 


1.60 


0.334 


S30WH07 


30 


13.8 


1.49 


0.219 


S35WH07 


35 


13.6 


1.61 


0.369 


S40WH07 


40 


15.3 


1.83 


0.599 


S45WH07 


45 


13.0 


1.79 


0.556 


S50WH07 


50 


9.76 


1.50 


0.221 


S60WH07 


60 


7.25 


1.46 


0.175 


S80WH07 


80 


6.33 


1.48 


0.210 


S100WH07 


100 


6.04 


1.46 


0.247 


S120WH07 


120 


5.96 


1.43 


0.172 


m35OCWH06 


35 


28.1 


2.08 


0.457 


E20HLW00 


20 


11.0 


1.74 


0.320 


E25HLW00 


25 


5.45 


1.70 


0.294 



Note. — The model name contains the information nec- 
essary to uniquely specify the presupernova model. For 
nonrotating progenitors, the beginning letter in the model 
name refers to th e metallicity of the pr ogenitor, following 
the convention of[Wooslev et al. (2002|), "s", "u", and "z" 
are used for solar, 10^* solar, and zero metalhcities, respec- 
tively. Following is the ZAMS mass, next we specify the 
progenitor model set (see the text for references). For ro- 
tating progenitors, we follow the naming convention of the 
original reference. 

a We define the iron core edge to be where Y e = 0.495. 

b €2 5 is determined at bounce in collapse runs using the 

LS180 EOS and will vary only slightly with EOS. 
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Figure 2. Presupernova mass as a function of ZAMS mass for the various 
model sets considered here. See the text for discussion. 



An additional uncertainty in massive s tar evolution is the 
phenomenon of large episodic mass loss dSmithll2007l) . Un- 
knowns and uncertainties in both the cause and effect of large 
episodic mass loss currently prevent detailed stellar evolution 
calculations from including this phenomenon at all. 

3.2. Grid Setup 

Based on resolution studies, we employ a computational 
grid setup with a total of 1050 zones. Near the origin and ex- 
tending out to 20 km, we employ a constant grid spacing of 
80m (250 zones). Outside of 20km we logarithmically space 
the remaining 800 zones to a radius where the initial density 
falls to 2000gcm~ 3 . We require the high resolution near the 
center for late postbounce times when the postshock region 
becomes small (r s h oc k ^ 20 km) and when the PNS is close to 
dynamical collapse to a BH. We interpolate the various pre- 
supernova profiles (p, T, Y e , v, fi) to our grid using linear 
interpolation. 

3.3. Rotation 

In simulations including 1.5D rotation, we directly use the 
angular velocity of the progenitor model if it was evolved with 
rotation or assign specific angular momentum via the rotation 
law 



j(r) = 7"i6,c 



1 + 



10 16 cm 2 s _1 , (11) 



where y'i6.oo is the specific angular momentum at infinity in 
units of 10 16 cm 2 s" 1 . We define Am to be the radius where 
the enclosed mass is 1M Q . This is a variation on the rotation 
law commonly used in simulations of rotating core collapse 
(e.g. JOtt et al.ll200"6l) . where O(r) = j(r)r~ 2 is prescribed and 
the differential-rotation parameter A is set to some constant 
radius. The advantage of prescribing j (which is conserved 
along Lagrangian trajectories) and choosing the value of A 
based on a mass coordinate is that progenitors from different 
groups that are evolved to different points still yield similar 
PNS angular momentum distributions for a given choice of 
ji6,oo- Equation ( fTTT ) leads to roughly uniform rotation in the 



core inside Am q (j(r) oc r 2 ) and angular velocity Q(r) decreas- 
ing with r 2 further out (J(r) = const). We note that when 1 M Q 
of material is contained within 10 3 km, which is typical of 
many progenitors, the central rotation rate is 7i6 j00 rads _1 . 

Our way of assigning rotation to precollapse models ap- 
proximates well the predictions of core rotati on (inner —few 
Mq) from stellar evolution studies (see, e.g., lOtt et alJl2006l 
for comparison plots) and, thus, is useful for studying rota- 
tional effects on BH formation. Equation ( TTTb does not, how- 
ever, capture the rise in specific angular momentum observed 
at larger radii (or mass coordinate) that is important for the 
potential evolution toward a l ong GRB and seen in rec ent ro- 
tating progenitor models (e.g. JWooslev & Hegerll200"6l) . 

4. RESULTS 

4.1. Fiducial Model 

We begin our discussion with a detailed description of the 
evolution of a failing CCSN from core collapse, through 
bounce, and the subsequent postbounce evolution to BH for- 
mation. For this, we choose the 40 M Q ZAMS-mass pro- 
genitor model s40WH07. We evolve this progenitor using 
the LSI 80 EOS, do not include rota tion, and use the stan- 
dard setting of /heat = 1 (see Section 12.2b . In Figure [3] we 
show the evolution of the radial coordinate of select baryonic 
mass shells as a function of time and we highlight shells en- 
closing 0.5, 1.0, 1.5, 2.0, and 2.5M Q . In addition, the fig- 
ure shows the shock radius and the positions of the energy- 
averaged v e and v x neutrinospheres as a function of time. 
The prebounce collapse phase (t < 0) lasts ~ 450ms. At 
bounce, the central value of the lapse function is a c ~ 0.82, 
and the metric function X has a maximum of — 1.1 and peaks 
off-center at a baryonic mass coordinate of ~ 0.6M Q which 
roughly corresponds to the edge of the inner core. The in- 
ner core initially overshoots to a maximum central density 
p c ~ 5.0 x 10 14 gcm" 3 , then settles at ~ 3.7 x 10 14 gcm" 3 . p c 
subsequently increases as accretion adds mass to the PNS. 
The bounce shock forms at a baryonic mass coordinate of 
<~ 0.6M Q . From there, it moves out quickly in mass, reaching 
a baryonic mass coordinate of ~ 1 .5 M at 22 ms after bounce, 
2Mq at ~ 162ms, and 2.25M at ~ 329ms. In radius, the 
shock reaches a maximum of ~ 120km at 38 ms after bounce. 
There it stalls, then slowly recedes. At 10 ms after bounce, the 
accretion rate through the shock is <~ 18M Q s" 1 and drops to 
-2.7,-1 .7, and - 1 .25 M Q s" 1 at 100, 200, and 300 ms after 
bounce, respectively. The drop in the accretion rate has lit- 
tle effect on the failing supernova engine. In agreement with 
previous w ork that employed a more accurate neutrino treat - 
ment (e.g., iThompson et al.l 120031: iLiebendorfer et all 120051) . 
the ID neutrino mechanism is manifestly ineffective in driv- 
ing the shock, yielding, in this model, a heating efficiency 
V = ^absorbed/i^+p,, of only — 3% (on average). The neutri- 
nospheres (where the energy-averaged optical depth t = 2/3) 
are initially exterior to the shock but are surpassed by the lat- 
ter in a matter of milliseconds after bounce, leading to the v e 
deleptonization burst. At all times, the v x neutrinosphere is 
interior to the v e neutrinosphere, which in turn is slightly in- 
terior to the v e neutrinosphere. The mean neutrino energies 
also follow this order. They are the largest for v x and the low- 
est f or v e and increase with decreasing neutrinosphere radii 
e.g. IThompson et alj|2003h ISumivoshi et al.ll2007fc lOtt et alj 
2008t jFischer et alj|2009h 
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Figure 3. Evolution of baryonic mass shells in the nonrotating model 
.S-40WH07 evolved with the LSI 80 EOS. We also include the shock location 
and the radii of the v e and v x neutrinospheres. The P c -sphere (not shown), 
is inside, but very close to the ^-sphere. The vertical dotted line denotes a 
change of timescale in the plot, highlighting the final ~ 1ms of evolution be- 
fore the central density reaches ~ 4.2 X 10 15 gem -3 and the simulation halts. 
We specifically highlight the 0.5, 1.0, 1.5, 2.0, and 2.5Mq baryonic mass 
shells with dashed lines. With solid lines, for M < 2Mq, we plot every 0.1 
Mq mass shell. Above 2Mq , we plot mass shells with a spacing of 0.05 Mq . 

At ~ 408 ms after bounce, the shock has receded to ~ 20 km 
and the PNS has reached a baryonic (gravitational) mass of 
~ 2.33M© (~ 2.23M Q ). The difference between baryonic 
and gravitational mass, at this point in the evolution, is due 
to the ~ 1.9 x 10 53 erg of energy radiated by neutrinos. At 
this point, dynamical PNS collapse to a BH sets in and hap- 
pens on a coordinate timescale of < 1 ms. In the rightmost 
part of Figure [3] we zoom in to the final 1 ms of evolution to 
show detail. The first signs of collapse manifest themselves 
in the development of a radial infall velocity profile at the 
PNS edge. The PNS then collapses in on itself and the central 
density increases by a factor of ~ 3 in only ~ 1 ms of coordi- 
nate time. The simulation crashes due to EOS limitations at 
p c - 4.2 x 10 15 gcm" 3 and with a c = 0.006. At this point the 
peak of the metric function X = [1— 2/n(r)/r] _1 / 2 is ~ 4.4 at a 
coordinate radius of ~ 6.8km. There, the fluid velocity also 
peaks at ~ -0.83c. The shock recedes by ~ 8km in the last 
~ 1 ms of evolution to a radial coordinate of ~ 12 km. During 
the last ~ 0.05 ms, due to the central lapse dropping to nearly 
zero, the evolution of the mass shells slows near the origin. 
This is characteristic for our choice of gauge. If the simula- 
tion were to continue, X would become singular at the event 
horizon that would appear after infinite coordinate time in our 
coordinates dPetrich et al.lll986l) . 

The .S-40WH07 model discussed here is a typical example 
of a failing CCSN in spherical symmetry. We present the re- 
sults of a large number of such models in Table |2j where for 
each EOS and progenitor model we show the time to BH for- 
mation as measured from bounce and the mass, both baryonic 
and gravitational, of the PNS when the central value of the 
lapse function a reaches 0.3 (roughly the point of instability). 
In this table, the model name describes the initial model. The 
metallicity is denoted by one of three letters: s, u, and z which 
represent solar, 10~ 4 solar, and zero metallicity, respectively. 
Following the metallicity is the ZAMS mass and the progen- 
itor model set. In many simulations, particularly in those 



employing stiff EOS, a BH does not form within 3.5s. For 
these simulations we include in parentheses the mass inside 
the shock at 3.5 s. We note that at BH formation the shock is 
typically at a distance of < 20 km and there is very little mass 
between the shock and the PNS. The dynamical collapse to a 
BH happens very quickly (t < 1 ms) during which very little 
additional accretion occurs. 

4.2. Comparison with Previous Work 

The ,?40WW95 progenitor was consider e d in the BH for- 
mation studies of lLieberid orfer et afj (120041) . ISumivoshi et al. 
( 20071) (hereinafter referred to as S07), and iFischer et al. 
(2009) (hereinafter referred to as F09). For comparison, we 
perform simulations with this progenitor for both the LSI 80 
and HShen EOS. Table [3] compares two key quantities, the 
time to BH formation and the maximum baryonic PNS mass, 
obtained with GRID with the results obtained in the afore- 
mentioned studies. 

For the LSI 80 EOS, we find a time to BH formation of 
~ 524ms and a maximum baryonic PNS mass of ~ 2.26M©, 
which is ~ 3 % larger than predicted by F09. We attribute 
this discrepancy to the different neutrino transport methods 
used. GRID'S leakage scheme has the tendency to somewhat 
over predi ct electron-type neutrin o luminosities (see the dis- 
cussion in |Q' Connor & Ottll2010l) . resulting in lower gravita- 
tional masses compared to full Boltzmann transport calcula- 
tions. Our time to BH formation is longer by ~ 100 ms or 
~ 20%. This disagreement is relatively larger than the bary- 
onic mass disagreement due to the low accretion rate at late 
times that translates small differences in mass to large differ- 
ences in time. At ~ 435.5ms, the time to BH formation of 
F09, our PNS has a baryonic mass of ~ 2.17 M©, which is 
consistent to ~ 1 % with F09. We find it more difficult to rec- 
oncile our results (and those of iLiebendorfer et"aT] (12004ft and 
F09) with the simulations of S07. Their maximum PNS bary- 
onic mass and the time to BH formation suggest a lower accre- 
tion rate throughout their evolution (~ 2.1 M© in ~ 560ms). 

In the simulation run with the stiffer HShen EOS, the larger 
maximum PNS mass leads to a delay of BH formation until a 
postbounce time ~ 1.129 s and we find a maximum baryonic 
PNS mass of ~ 2.82M Q . The maximum PNS mass and time 
to BH formation of S07 again suggest an accretion rate in dis- 
agreement with F09 and our work. The results of F09 with 
the HShen EOS suffer from a glitch in F09's EOS table inter- 
polation scheme which has since been fixed (T. Fischer 2010, 
private communication). This leads to a postbounce time to 
BH formation of ~ 1 .4 s and a maximum baryonic PNS mass 
of ~ 3.2M Q . Results from more recent simulations correct 
this error and are presented in Table [3] (T. Fischer 2010, pri- 
vate communication). 

4.3. Equation- of- state Dependence and Thermal Effects 

The maximum PNS mass and, thus, the evolution toward 
BH formation, d epends strongly on the EOS. This was re- 
alized early on dBurrowsl 119881) and has recently been in- 
vestigated by S07 and F09 who compared models evolved 
with the LSI 80 and HShen EOS. Here we extend their dis- 
cussion and include also the LS220 and LS375 EOS. For a 
given accretion history, set by progenitor structure and in- 
dependent of the high-density EOS, a stiffer nuclear EOS 
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iw b.max 


M 
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Ml 
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M 
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.S-20WW95 


0.787 


2.238 


2.108 


1.129 


2.377 


2.201 


3.351 


3.060 


2.653 


2.287 


2.751 


2.486 


.«25WW95 


0.737 


2.246 


2.118 


1.046 


2.383 


2.211 


2.707 


3.054 


2.656 


1.990 


2.760 


2.498 


s40WW95 


0.524 


2.263 


2.137 


0.666 


2.406 


2.240 


1.381 


3.043 


2.674 


1.129 


2.815 


2.562 


.S20WHW02 




(1.949) 


(1.794) 




(1.950) 


(1.798) 




(1.951) 


(1.807) 




(1.943) 


(1.805) 


.V25WHW02 


1.021 


2.211 


2.079 


1.504 


2.355 


2.172 




(2.917) 


(2.559) 


2.929 


2.736 


2.468 


.S30WHW02 


1.820 


2.144 


1.978 


2.986 


2.331 


2.108 




(2.416) 


(2.182) 




(2.405) 


(2.190) 


.«35WHW()2 


2.073 


2.141 


1.976 


3.334 


2.328 


2.105 




(2.351) 


(2.137) 




(2.340) 


(2.141) 


.V40WHW02 


1.512 


2.168 


2.019 


2.231 


2.336 


2.134 




(2.634) 


(2.355) 




(2.615) 


(2.362) 


.«75WHW02 




(1.920) 


(1.781) 




(1.920) 


(1.784) 




(1.921) 


(1.791) 




(1.913) 


(1.787) 


1/20WHW02 


0.938 


2.215 


2.082 


1.367 


2.358 


2.175 




(2.852) 


(2.516) 


3.004 


2.734 


2.466 


U25WHW02 


1.759 


2.160 


2.009 


2.798 


2.330 


2.124 




(2.446) 


(2.218) 




(2.429) 


(2.220) 


»30WHW02 


0.922 


2.217 


2.084 


1.353 


2.359 


2.178 




(2.802) 


(2.483) 


3.228 


2.731 


2.462 


»35WHW02 


0.379 


2.347 


2.242 


0.484 


2.465 


2.329 


1.308 


3.020 


2.693 


1.075 


2.847 


2.608 


U40WHW02 


0.369 


2.346 


2.241 


0.453 


2.469 


2.333 


0.946 


3.023 


2.710 


0.849 


2.874 


2.638 


U45WHW02 


0.441 


2.301 


2.187 


0.548 


2.433 


2.284 


1.108 


3.027 


2.694 


0.959 


2.842 


2.600 


1/50WHW02 


0.563 


2.273 


2.154 


0.706 


2.408 


2.251 


1.365 


3.030 


2.676 


1.163 


2.816 


2.569 


U60WHW02 


0.432 


2.363 


2.267 


0.579 


2.460 


2.331 


1.346 


3.009 


2.693 


1.165 


2.849 


2.620 


M75WHW02 


0.226 


2.526 


2.449 


0.285 


2.592 


2.498 


0.626 


3.006 


2.775 


0.594 


2.984 


2.791 


Z20WHW02 


3.295 


2.116 


1.934 




(2.141) 


(1.955) 




(2.143) 


(1.968) 




(2.132) 


(1.968) 


Z25WHW02 


0.602 


2.283 


2.167 


0.956 


2.398 


2.239 


3.443 


3.050 


2.650 


2.351 


2.762 


2.505 


Z30WHW02 


1.772 


2.149 


1.989 


2.964 


2.329 


2.114 




(2.413) 


(2.187) 




(2.401) 


(2.192) 


Z35WHW02 


0.446 


2.321 


2.213 


0.619 


2.436 


2.291 


1.939 


3.027 


2.669 


1.380 


2.813 


2.569 


Z40WHW02 


0.365 


2.350 


2.245 


0.450 


2.471 


2.335 


0.958 


3.023 


2.711 


0.856 


2.874 


2.639 


.«25LC06A 


1.220 


2.176 


2.029 


2.547 


2.333 


2.130 




(2.440) 


(2.213) 




(2.398) 


(2.195) 


.V30LC06A 


1.101 


2.181 


2.035 


1.726 


2.342 


2.141 




(2.767) 


(2.446) 




(2.695) 


(2.421) 


.V35LC06A 


1.029 


2.186 


2.040 


1.726 


2.338 


2.133 




(2.567) 


(2.305) 




(2.517) 


(2.285) 


.V40LC06A 


0.746 


2.232 


2.102 


1.138 


2.372 


2.193 




(2.796) 


(2.470) 


3.390 


2.723 


2.452 


,v60LC06A 


0.393 


2.331 


2.224 


0.512 


2.450 


2.310 


1.536 


3.025 


2.678 


1.278 


2.816 


2.572 


.S80LC06A 


0.429 


2.308 


2.197 


0.530 


2.437 


2.293 


1.075 


3.021 


2.689 


1.083 


2.825 


2.581 


5120LC06A 


0.262 


2.439 


2.351 


0.317 


2.531 


2.423 


0.661 


3.001 


2.745 


0.728 


2.911 


2.701 


440LC06B 


0.958 


2.189 


2.043 


1.411 


2.349 


2.152 




(2.957) 


(2.576) 


2.887 


2.720 


2.444 


.S-60LC06B 


3.073 


2.117 


1.934 




(2.166) 


(1.972) 




(2.165) 


(1.984) 




(2.126) 


(1.961) 


.S-80LC06B 


2.441 


2.131 


1.963 




(2.260) 


(2.052) 




(2.264) 


(2.071) 




(2.249) 


(2.069) 


sl20LC06B 


2.983 


2.120 


1.944 




(2.171) 


(1.984) 




(2.167) 


(1.992) 




(2.102) 


(1.947) 


.V20WH07 


1.275 


2.180 


2.035 


1.876 


2.341 


2.143 




(2.712) 


(2.412) 




(2.694) 


(2.426) 


s25WH07 


1.066 


2.202 


2.065 


1.523 


2.352 


2.165 




(2.975) 


(2.595) 


2.796 


2.736 


2.466 


.V30WH07 


1.751 


2.150 


1.991 


2.978 


2.329 


2.115 




(2.408) 


(2.184) 




(2.397) 


(2.190) 


.V35WH07 


0.836 


2.232 


2.104 


1.203 


2.369 


2.194 




(2.918) 


(2.563) 


2.689 


2.744 


2.481 


s40WH07 


0.408 


2.334 


2.228 


0.561 


2.448 


2.306 


1.596 


3.024 


2.680 


1.259 


2.827 


2.585 


s45WH07 


0.454 


2.319 


2.210 


0.626 


2.435 


2.289 


2.027 


3.028 


2.667 


1.395 


2.812 


2.567 


.V50WH07 


1.813 


2.147 


1.987 


2.989 


2.329 


2.113 




(2.411) 


(2.185) 




(2.399) 


(2.190) 


s60WH07 


2.778 


2.124 


1.947 




(2.230) 


(2.023) 




(2.232) 


(2.039) 




(2.220) 


(2.040) 


.V80WH07 


2.113 


2.139 


1.974 


3.284 


2.328 


2.104 




(2.363) 


(2.145) 




(2.350) 


(2.148) 


4-100WH07 


1.457 


2.163 


2.008 


2.355 


2.335 


2.124 




(2.539) 


(2.281) 




(2.524) 


(2.289) 


i'120WH07 


3.043 


2.120 


1.940 




(2.179) 


(1.985) 




(2.180) 


(1.999) 




(2.169) 


(1.999) 



Note. — BH formation times and maximum PNS mass (both baryonic and gravitational) for nonrotating runs with /beat = 1 for all four EOS. 
We stop our simulations at 3.5 s after core bounce. Models that have not formed a BH by then probably explode in nature. They are marked by 
• • • , but we include the masses insi de the shoc k at 3.5 s in pare ntheses. The progenitor models are the r esult of vario us stellar evolution studie s: 
WW95, Woosley & Weaye^ fT995T) ; WHW02. IWooslev et alj J2003) ; LC06. ILimongi"& Chieffi 1 200|); and WH07, |Woosley & Hegerj <2007T) . 
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Table 3 

A-40WW95 Comparison with LSI 80 and HShen EOS. 



LSI 80 



HShen 



Liebendorfer et al. (2004) 
Sumivoshi et al. (2007) 
Fischer et al. (2009) 
This work 









M fc>max 




(s) 


(M @ ) 


o) 


(Mq) 


's 


~0.5 


-2.20 






U 


0.56 


2.1 


1.34 


2.66 


bfj 


0.4355 


2.196 


1.030 a 


2.866 a 




0.524 


2.263 


1.129 


2.815 


O 
< — i 



a See the text fo r a discussion of the HShen EOS results from 
IFischer et aTH2009l) . 

leads to a larger postbounce time to BH formation. In Fig- 
ure |U we plot the evolution of the central density p c of the 
s40WH07 model evolved with the four considered EOS. Each 
EOS leads to a characteristic maximum central density at 
bounce that is practically independent of progenitor model: 

3.7 x 10 14 gcm"\ 
and 

HShen EOS, respectively. As expected, the variant using the 
softest nuclear EOS (LSI 80) shows the steepest postbounce 
increase in p c and becomes unstable to BH formation at only 
<~ 408 ms for this progenitor. The onset of BH formation is 
marked by a quick rise in the central density. This is most ob- 
vious from the p c evolutions of the model variants using the 
stiff HShen and LS375 EOS. 



-4.8 x 10 14 gcm" 3 , 



4.4x 10 14 gcnT 3 , 



and - 3.4 x 10 14 gcm" 3 for the LS180, LS220, LS375 



Interestingly, the nominally stiffest EOS (LS375) leads to 
higher central densities than the softer HShen EOS up un- 
til ~ 1.1s after bounce. We find that this is due to the 
HShen EOS yielding higher pressure at p < 3 x 10 14 gcm~ 3 , 
T ~ lOMeV, and Y e ~ 0.3. This higher pressure, initially in 
the core and later in the outer PNS layers, maintains the PNS 
at a lower central density. The cold-NS mass-radius relation 
shown in Figure[T]also exhibits this. For a given low-mass NS, 
the HShen EOS predicts a lower central density. For cold NSs, 
this trend continues until p c ~ 5.4 x 10 14 gcm" 3 . Thermal ef- 
fects, which are stronger in the HShen EOS, will increase this 
value for hot PNSs. 

We also plot in Figure|4]the evolution of the mass accretion 
rate M in model .s40WH07 (evaluated at a radius of 200 km). 
Variations in the high-density EOS have no effect on M which 
is most sensitive to progenitor structure. Sudden drops in M 
occur when density discontinuities that go along with com- 
positional interfaces advect in. An example of this can be 
seen in s40WH07 at ~ 400 ms after bounce where M drops 
by ~ 30 % due to a density jump at a baryonic mass coor- 
dinate of ~ 2.35M Q . Su ch interfaces are com mon features 
of evolved massive stars dWooslev et al .1120021) and can help 
jumpstar t shock revival in spe cial cases (see , e.g ., the 1 1.2M Q 
model of iBuras et al.l d2006al) . and Section 1431 of this study). 
In the BH-formation context, they lead to a disproportionate 
increase in the time to BH formation in models whose EOS 
permit a PNS with mass greater than the mass coordinate of 
the density jump. 

The maximum gravitational (baryonic) PNS masses for the 
four models shown in Figure|4]are ~ 2.23M Q (~ 2.33M Q ), 
- 2.31M Q (~ 2.45M©), - 2.68M Q (~ 3.O2M ), and ~ 
2.59M Q (— 2.83M ); and the BH formation times are ~ 
408ms, - 561ms, - 1.596s, and - 1.259s for the LSI 80, 
LS220, LS375, and HShen EOS, respectively. The maximum 
cold NS gravitational masses are ~ 1.83M Q , ~ 2.04 M Q , 
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, M, 
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Figure 4. Central density (left ordinate) and accretion rate (right ordinate) vs. 
time since bounce for the .s40WH07 progenitor and four EOS. BH formation 
occurs when the central density diverges. Each p c curve is annotated with the 
maximum gravitational PNS mass. The drop in the accretion rate at ; ~ 0.4 s 
is due to the accretion of a mass shell where the density drops by ~ 30%. 
Note the accretion rate is on a logarithmic scale. 



~ 2.72M , and - 2.24M Q for the LS180, LS220, LS375, 
and HShen EOS, respectively. 

In models evolved with the LSI 80, LS220, and HShen 
EOS, the maximum gravitational PNS mass is larger than the 
maximum gravitational cold NS mass. We can understand the 
differences between these cold NS and PNS maximum masses 
by comparing the PNS structure with various TOV solutions. 
In Figure [5] we plot the density and temperature profiles of 
the s40WH07 model evolved with the HShen EOS just prior 
to collapse to a BH. At this time, t ~ 1 .098 s, the central lapse 
is a c = 0.35, the central density is p c ~ 1.44 x 10 15 gcm" 3 , 
T c ~ 42. 4 Me V, and the PNS gravitational (baryonic) mass 
is ~ 2.5 1M Q (~ 2.74M Q ). For comparison, we include in 
Figure [3] three TOV solutions, all with the same central den- 
sity but different temperature and Y e profiles. Specifically, we 
plot the density profile assuming (1) T(r) = O.lMeV, this is 
the "cold" NS case, (2) T(r) = 42.4 Me V, which corresponds 
to the central temperature from the GRID evolution, and (3) 
T(r) = r GR i D , assuming the same radial temperature profile as 
the GRID model. We impose neutrinoless (5 equilibrium for 
the former two TOV solutions and, similar to the temperature, 
assume the Y e profile of the GRID model for the latter. For 
this comparison, GRID is run without neutrino pressure and 
energy contributions, since they are also neglected in the TOV 
solution. 

Inside of ~ 6 km, corresponding to a gravitational mass co- 
ordinate of ~ 0.4Mq, the material is not shock heated but 
rather is heated only via adiabatic compression. The outer 
regions (^6-11 km) of the PNS are hot compared to the in- 
ner core. This is due to accretion and compression of shock 
heated material onto the PNS surface. In this region, the ther- 
mal pressure support is sufficiently strong to flatten out the 
PNS density profile. More mass is located at larger radii 
compared to constant-temperature TOV solutions. This de- 
creases PNS compactness, increasing the maximum gravita- 
tional mass. The cold-NS and the T = T C TOV solutions have a 
gravitational mass of ~ 2.23M Q and ~ 2.35M Q , respectively. 
On the other hand, the TOV solution that assumes the same T 
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and Y e profile as the GRID model yields a gravitational mass 
of ~ 2.46M , within ~ 2% of the PNS gravitational mass 
in the full GRID simulation. Tests in which we vary the Y e 
distribution in the TOV solution with T = T(r) show that the 
maximum PNS mass is insensitive to variations in Y e from 
the GRID profile to neutrinoless /^-equilibrium. All this leads 
us to the conclusion that it is thermal pressure support in the 
outer PNS core that is responsible for increasing the maxi- 
mum stable gravitational PNS mass beyond that of a cold NS. 
Our findin g is in agreement with the rece nt BH formation sim - 
ulations of ISumiyoshi et all (120071) and iFischer et alJ (12009ft . 
who noted the same differences to cold TOV solutions, but 
did not pinpoint th eir precise cause . However, our result is in 
disagreement with iBurrowsl dl988l) who reported maximum 
PNS masses within a few percent of a solar mass off their 
cold-NS values. This could be related to Burrows's specific 
choice of EOS. As we discuss below, stiff nuclear EOS have a 
more limited response to thermal effects. Another resolution 
to this disagreement could be the nature of his PNS cooling 
simulations that were not hydrodynamic, but rather employed 
a Henyey relaxation approach with imposed accretion. 

We find the same overall systematics of increased maxi- 
mum PNS mass due to thermal pressure support for the en- 
tire set of progenitors evolved with the LSI 80, LS220, and 
HShen EOS (variations due to diffe rences in progenitor struc- 
ture are discussed in Section 14. A\ . In the sequence of the 
LS EOS, the relevance of thermal pressure support decreases 
with increasing stiffness. In the case of the perhaps unphysi- 
cally stiff LS375 EOS, the effect of high temperatures in the 
outer PNS core is reversed: the PNSs in GRID simulations 
become unstable to collapse at lower maximum masses than 
their cold counterparts. This very surprising observation is 
understood by considering that in GR, higher temperatures 
not only add thermal pressure support to the PNS, but also in- 
crease its mass-energy. This results in a deeper effective po- 
tential well and, thus, is destabilizing. In the LSI 80, LS220, 
and HShen case, the added thermal pressure support is sig- 
nificant and dominates over the latter effect. In the very stiff 
LS375 EOS, the added thermal pressure component is negli- 
gible, and the destabilizing effect dominates. 

Finally, we point out quantitative differences in models 
evolved with and without neutrino pressure in the dense 
neutrino-opaque core. In the s40WH07 model evolved with 
the HShen EOS, the difference in the maximum gravitational 
mass is <~ 0.08 M Q (~ 3 %) and the difference in the time to 
BH formation is ~ 160ms (<~ 14%). These numbers depend 
on the employed EOS and progenitor model. In test calcula- 
tions with a variety of progenitors and EOS, we generally find 
increases of the maximum PNS gravitational mass of < 5 %. 

4.4. Influence of Pre supernova Structure 

The failure of a CCSN becomes definite only when accre- 
tion pushes the PNS over its maximum mass and a BH forms. 
Hence, the time to BH formation is a hard upper limit to the 
time available for the supernova mechanism to reenergize the 
shock. We will demonstrate in this section that it is possible 
to estimate, for a given nuclear EOS, the postbounce time to 
BH formation in non- or slowly spinning on the basis of a 
single parameter, progenitor bounce compactness, £2.5, which 
we introduced in Section [3] In Figure [6] we plot the post- 
bounce time to BH formation (Teh) as a function of £2.5 for 
all nonrotating models considered in this study and listed in 
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Figure 5. Comparison of radial density (left ordinate) and temperature (right 
ordinate) profiles of the PNS just before collapse (a c = 0.35) in model 
.S-40WH07 evolved with the HShen EOS with profiles obtained from a TOV 
solution using the same central density and the same radial temperature and 
Y e distributions as in model J40WH07 (dashed). For comparison, we also 
include profiles obtained with the TOV equations assuming both T = T c = 
42.4 MeV (dot-dashed) and T = 0.1 MeV "cold" (dot-dot-dashed) and /3- 
equilibrium. The flattening of the density profile between 5 and 1 1 km is 
due to the strong thermal pressure support in this region (dotted). The gravi- 
tational mass inside the shock (whose position we denote with a vertical black 
line) of the s40WH07 model and of the TOV star agree to within 2%. For 
this comparison, we switched off neutrino contributions to the internal energy 
and pressure in GRID. 



Table 12 The distribution of data points for each EOS can be 
fit with a function oc (£2.5) 3//2 - This remarkable result can be 
understood as follows: using Kepler's third law, consider the 
Newtonian free fall time to the origin for a mass element dm 
initially located at r* and on a radial orbit about a point mass 
of M* > dm, 



t& '■ 



GM* 



rl 



8GM* 



(12) 



Here, for clarity, the quantities are in cgs units. G is the 
gravitational constant and a is the semimajor axis equal to 
half of the apoapsis, r*. Recalling the definition of £2.5, if the 
mass element dm is located at a mass coordinate of 2.5 M© 
and at a radial coordinate of r„, then r„ = 2500 km/^. 5 , and 
we can write the free fall time in terms of £2.5, 



.2.5M, 



£-"" =O.241(&5) 



-3/2 



(13) 



In Figure[6] we overplot this Newtonian free fall time for a 
mass element at baryonic mass coordinate 2.5 Mq as a func- 
tion of £2.5- For small £2.5, the mass element begins its free 
fall from a large radius and, hence, takes longer to reach to 
origin. In general, material in outer layers of the star will not 
begin to fall freely until it loses pressure support. Hence, the 
free fall appro ximation is not exact (within a factor of ~ 2; 
IBurrowsl 19861) . but describes the general behavior of ?bh very 
well. The deviation of data points from the free fall curve is 
because the maximum PNS mass is different for each model 
and EOS. Models evolved with the LS180 EOS have PNSs 
with maximum baryonic masses ranging from 2.1 to 2.5M Q . 
Models with low £2.5 correspond to the lower end of this mass 
range. For these models, ?bh can be somewhat less than the 
free fall time of the 2.5 Mq mass element, because less mate- 
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Figure 6. BH formation time as a function of the bounce compactness (£2.5) 
for all nonrotating models presented in Table|2]that form BHs within 3.5 s of 
bounce. Simulations performed with the LSI 80, LS220, LS375, and HShen 
EOS are labeled with circles, squares, diamonds, and triangles, resp ecti vely. 
Also shown (dashed line) is the free fall time to the origin (Equation ll3t of a 
mass element located at a baryonic mass coordinate of 2.5 Mq. 



rial is needed to form a BH. The maximum baryonic PNS 
mass range for models using the LS220 EOS is somewhat 
higher, 2.3-2.6M©. BH formation times for these models 
are more in line with the Newtonian free fall prediction. Mod- 
els evolved with the LS375 and HShen EOS have PNSs that 
must accrete upward of ~ 3M Q of material before becoming 
unstable. This significantly increases ?bh above the free fall 
prediction for the £2.5 characteristic mass element. 

Thermal pressure support can increase the maximum grav- 
itational PNS mass (M g max ) as we have seen in Section [431 
for the s40WH07 model. In Figure [7] we plot M g max as a 
function of £2.5 for all nonrotating models listed in Table [2] 
As obvious from this figure, M gmax depends in a predictable 
way not only on the EOS, but also on the bounce compact- 
ness of the presupernova model. Progenitors with high £2.5, in 
addition to forming BHs faster, create PNSs that are stable to 
higher masses. This is a simple consequence of the fact that 
progenitors with larger £2.5 have iron cores with systemati- 
cally higher entropies and masses signi ficantly above the cold 
Chandrase khar mass (see Table [T] and iBaron & Coopersteinl 
[19901 and E 

urrows & Lattimerl 119831) . Adiabatic collapse 
leads to higher PNS temperatures after bounce in progenitors 
with high £2.5 and, hence, more thermal support. This leads to 
higher maximum PNS masses. This effect can be large, up to 
25% for models with large £2.5 and soft EOS. 

4.5. Preventing BH Formation with Artificial 
Neutrino-driven Explosions 

In a successful CCSN, the shock is reenergized before 
enough material can accrete onto the PNS to make it unstable. 
While fully self-consistent spherically symmetric simulations 
generally fail to explode in all but a few very low mas s pro- 
genitors (cf. iKitaura et al. 20(3 iBurrows et al.ll2007ah . one 
can explode any star by the ID neutrino mechanism by arti- 
ficially increasing the energy deposition in the postshock re- 
gion. Without such an increase, all of our simulations fail to 
explode. Our parameterized heating (/heat in Equation |9]l al- 
lows us to explore "how much" neutrino heating is needed to 
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Figure 7. Maximum gravitational PNS masses as a function of the bounce 
compactness (§2.5) for all nonrotating models presented in Table[2]that form 
BHs within 3.5 s after bounce. Simulations performed with the LS180, 
LS220, LS375, and HShen EOS are labeled with circles, squares, diamonds, 
and triangles, respectively. Also shown (dotted lines, labeled) are the maxi- 
mum gravitational cold-neutron star (CNS) masses, M™ s ax , numerical values 
are 1.83, 2.04, 2.72, and 2.24M for the LS180, LS220, LS375, and HShen 
EOS, respectively. 



explode a given model (in ID). By comparison with results 
from previous self-consistent radiation hydrodynamics simu- 
lations we can then estimate whether a given progenitor and 
EOS combination is more likely to lead to an explosion or BH 
formation. 

Our method fo r drivi ng explosions is similar to 
iMurphy & Burrowsl (120081) . but has the advantage of 
being proportional to the neutrino luminosity obtained from 
the neutrino leakage scheme and therefore conserves energy. 
We iteratively determine the critical value of /heat to within 
1 % to what is needed for a successful explosion for a large 
subset of our models and the LSI 80, LS220, and HShen EOS. 
Of particular interest in this analysis is the time-averaged 
heating efficiency of the critical model, ?7hear ^ e define ^h eat 
as 



?7heat — 



41 dV 



(L„„+Lc 



(14) 



where ql is the net energy deposition rate and the neutrino 
luminosities are taken at the gain radius. We perform the time 
average between bounce and explosion, the latter time defined 
as when the postshock region assumes positive velocities and 
accretion onto the PNS ceases. is a useful quantity be- 
cause it characterizes how much of the available luminosity 
must be redeposited on average to explode a given progenitor. 
This is rather independent of transport scheme and code. For 
example, for the 15 Mq ZAMS solar- me tallicity progenitor of 
IWooslev & Weaver! ( fl995l) we find 7^* - 0.13. iBuras et alj 
d2006bl) who also artificially exploded this progenitor in ID, 
though with much more sophisticated neutrino transport, fincQ 
an average heating efficiency of 0. 1- 0.15 which is consisten t 
with our result. Note, however, that iMarek & Jankal (120091) 
observed in the same progenitor the onset of a self-consistent 
neutrino-driven explosion in 2D at an average heating effi- 



2 This we deduce from their Figure 28, bottom panel. 
<5,£ coo l includes all neutrinos, not just u e and v e . 



Note that their 
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Figure 8. r^"^ obtained with GRID as a function of bounce compactness. 
Plotted are models from the sWH07 data set using the LSI 80, LS220, and 
HShen EOS; and models from the wWHW02 data set using the LS220 EOS. 



ciency of ~ 0.07. This indicates a dependence of on 
dimensionality, should be kept in mind, and is consistent with 
recent work that suggest that dimension ality may be the key 
to successful neutrino-drive n explosions dMurphy & Burrows! 
l2008tlNordhaus et alJl2010h . 

Since GRID'S leakage/heating scheme is only a rough ap- 
proximation to true neutrino transport, and because our simu- 
lations assume spherical symmetry, we cannot make very ro- 
bust quantitative predictions for any one particular model, but 
rather study the collective trends exhibited by the entire set of 
62 progenitors that we consider here. In Figure [8] as a func- 
tion of bounce compactness £2.5, we plot fj^ t for all consid- 
ered models and EOS. The data are summarized in Table H] 
We can divide the results into two general regimes: models 
with £2.5 < 0.45 and those with £2.5 > 0.45. 

For many models with £2.5 ^ 0.45, oscillations in the shock 
position are ubiquitous near the transition from failing to 
exploding supernovae in ID (cf. iMurphv & Burrows! 120081: 
iBuras et al.ll2006bl: iFernandez & Thompsonll2009l) ." For both 
the LSI 80 and LS220 EOS, the jjg* required for an explo- 
sion, modulo noise, is roughly constant and ~ 0. 16 on average 
for low £2.5 models. Hence, explosion is the likely outcome 
of core collapse for progenitors with £2.5 ^ 0.45 if the nuclear 
EOS is similar to the LSI 80 or LS220 case. 



Table 4 

Explosion Properties 



The noise in the reheat distribution (absolute variations by 
up to ~ 10%) is in part a consequence of variations in post- 
bounce dynamics, such as the number and duration of pre- 
explosion oscillations. Compositional interfaces in some pro- 
genitor models, where jumps in the density lead to jumps in 
the accretion rate, also affect individual models leading to 
variations in fg* . For the LS180 and LS220 EOS, any differ- 
ences in with choice of EOS are indistinguishable given 
the noise in the data. 



For progenitors with £2.5 ^ 0.45, the f)h eat required to cause 
an explosion increases with £2.5 when run with the LSI 80 or 
LS220 EOS. Progenitors in this regime have tremendous post- 
bounce accretion rates, accumulating > 2M Q of baryonic ma- 
terial behind the shock within the first ~ 200ms after bounce. 
Without explosion, they form BHs within < 0.8 s (with the 
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0.409 


1.16 


0.163 


18 


0.195 


1.17 


0.188 


25 


0.334 


1.13 


0.158 


19 


0.177 


1.24 


0.146 


27 


0.258 


1.18 


0.153 


20 


0.288 


1.15 


0.176 


30 


0.219 


1.16 


0.179 


21 


0.143 


1.34 


0.133 


35 


0.369 


1.14 


0.164 


22 


0.292 


1.15 


0.181 


40 


0.599 


1.32 


0.266 


23 


0.453 


1.17 


0.165 


45 


0.556 


1.26 


0.245 


24 


0.410 


1.15 


0.163 


50 


0.221 


1.18 


0.172 


25 


0.334 


1.14 


0.185 


80 


0.210 


1.22 


0.142 


26 


0.235 


1.21 


0.142 




.S-WH07 HShen 




27 
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1.20 


0.152 


15 
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1.30 


0.245 


28 


0.274 


1.16 
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21 
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1.50 
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29 
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0.138 


23 
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1.27 


0.245 


30 


0.219 


1.18 
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24 


0.406 


1.31 


0.245 


31 


0.219 


1.21 


0.144 


25 


0.333 


1.49 


0.217 


32 


0.255 


1.17 


0.166 


27 


0.258 


1.52 


0.186 


33 


0.287 


1.15 


0.162 


30 


0.218 


1.32 


0.194 


35 


0.369 


1.13 


0.166 


35 


0.367 


1.37 


0.167 


40 


0.600 


1.30 


0.259 


40 


0.581 


1.22 


0.245 


45 


0.557 


1.25 


0.228 


45 


0.542 


1.24 


0.240 


50 


0.221 


1.19 


0.170 


50 


0.221 


1.41 


0.218 


55 


0.238 


1.24 


0.129 


80 


0.210 


1.50 


0.226 


60 


0.175 


1.29 


0.142 










70 


0.234 


1.21 


0.161 




sWW95 


LSI 80 




80 


0.210 


1.24 


0.143 


15 


0.088 


1.33 


0.130 


LOO 


0.247 


1.15 


0.175 










120 


0.172 


1.25 


0.152 
















HWHW02 LS220 








20 


0.338 


1.13 


0.155 


40 


0.721 


1.44 


0.297 


25 


0.223 


1.16 


0.168 


45 


0.656 


1.22 


0.267 


30 


0.326 


1.13 


0.156 


50 


0.575 


1.09 


0.174 


35 


0.666 


1.37 


0.284 


60 


0.624 


1.12 


0.133 


Note. 


— ^hc'i't colTes P on ds to the critical value needed to 


cause a 



successful explosion in GRID, jjj^* is the associated critical average 
heating efficiency defined in Equation U4t . 

LSI 80 and LS220 EOS). Hence, a very high heating effi- 
ciency of 0.23-0.27 is necessary to drive an explosion 
at early times against the huge ram pressure of accretion. It 
appears unlikely, even when multi-dimensional dynamics are 
factored in, that progenitors with £2.5 ^ 0.45 can be exploded 
via the neutrino mechanism. The most likely outcome of core 
collapse in such stars is BH formation. 

We draw the reader's attention to two outliers in the 
mWHW02 data set included in Figured the t/50WHW02 and 
m60WHW02 progenitors. These models have high £2.5, but 
feature compositional interfaces where the density drops by 
~ 50%. These are located at a mass coordinate of 1.82M© 
and 2.22 M Q in m50WHW02 and w60WHW02, respectively. 
When such an interface advects through the shock, the accre- 
tion rate drops suddenly, but the core neutrino luminosity re- 
mains large and an explosion is immediately launched. This 
results in a small value of f^ t and, therefore, in a low re- 
quired f/h eat . This demonstrates that the single parameter £2.5 
is not always sufficient to predict a progenitor's fate. 

In models with £2.5 ^ 0.45 and calculated using the HShen 
EOS, both fjfel t and are systematically higher than with 
the LSI 80 and LS220 EOS and explosion is less likely. Fur- 
thermore, the qualitative behavior of our simulations is dif- 
ferent with the HShen EOS. In many models with subcritical 
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Aeat and reheat, the shock is revived and begins to propagate 
to large radii of 0(1O 3 - 10 4 km), but the material behind it 
fails to achieve positive velocities. Hence, accretion onto the 
PNS is slowed but does not cease. High values of /heat are 
needed to avoid this and achieve full explosions. We caution 
the reader that this regime may not be well modeled by our 
neutrino treatment. Nevertheless, our results suggest that sys- 
tematically higher and fj"^ t are required to explode mod- 
els with the HShen EOS, even at low £2.5- In contrast to mod- 
els using the LSI 80 or LS220 EOS, models with £2.5 >0.45 
with the HShen EOS require roughly constant i^eat to explode. 
Since the HShen EOS can support a high maximum mass, the 
PNS can withstand BH formation longer and explosions may 
set in at later postbounce times when the accretion rate has 
dropped sufficiently. 

Finally, as an interesting aside, we point out the evolution 
of the M75WHW02 progenitor evolved with the LS220 EOS. 
This model has a bounce compactness of ~ 1.15 and, in the 
absence of an explosion, forms a BH ~ 0.285 s after bounce 
(with the LS220 EOS). This progenitor has a compositional 
interface at which the density drops by ~ 50% that is lo- 
cated at a baryonic mass coordinate of ~ 2.5 M . This is 
very close to the maximum mass of the m75WHW02 PNS 
(with the LS220 EOS), and well above the maximum cold 
NS (baryonic) mass. The model can be made to explode with 
/heat = 1-35 with a corresponding = 0.287. The result- 
ing PNS has a baryonic (gravitational) mass of ~ 2.54M© 
(2.44Mq). Interestingly, within ~ 100ms after the launch of 
the explosion, cooling of the outer PNS layers removes suf- 
ficient thermal pressure, rendering the PNS unstable to col- 
lapse and BH formation. This scenario will necessarily occur 
within the cooling phase for any PNS that is initially thermally 
supported above the maximum cold NS baryonic mass and is 
another avenue to BH formation. In our simulations, this con- 
dition is also met only in very few other models with very 
high £2.5 and fairly soft EOS, such as the 23, 40, and 45 M© 
progenitors from the sWH07 series using the LSI 80 EOS. In 
order to fully investigate this BH formation channel, a more 
sophisticated neutrino treatment is require d that allows accu - 
rate long-term modeling of PNS cooling dPons et al.lll999l) . 
since, in general, the Kelvin-Helmholtz cooling phase of PNS 
is O(10- 100s). 

4.6. Connection to Stellar Evolution and ZAMS Conditions 

4.6. 1 . ZAMS Mass and Metallicity 

Having established the systematic dependence of core col- 
lapse and BH formation on progenitor bounce compactness 
in Section 14.41 we now go further and attempt to connect to 
the conditions at ZAMS. Doing this is difficult, and, given 
the current state and limitations of stellar evolution theory 
and modeling, can be done only approximately. In general, 
presupernova stellar structure will depend not only on initial 
conditions at ZAMS (mass, metallicity, rotation), but also on 
particular evolution history and physics (binary effects, [ro- 
tational] mixin g, magnetic fields, nu clear reaction rates, and 
mass loss; cf. IWoosley et al.l 120021) . While keeping this in 
mind, we limit ourselves in the following to the exploration of 
single-star, nonrotating progenitors without magnetic fields. 
We focus on collapse models run with the LS220 EOS, but 
the general trends with EOS observed in the previous sections 
extend to here. 
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Figure 9. Bounce compactness (top panel) and time to BH formation (bot- 
tom panel) as functions of ZAMS mass for various p rog enitor sets. £2.5 is de- 
termined for each model at bounce using Equation 4 1 i . fBH for each model 
is obtained using the LS220 EOS and assuming no explosion. The times to 
BH formation for progenitor models that take longer than 3.5 s to form a BH 
are not shown. Breaks in the lines connecting models indicate this. For clar- 
ity, the time to BH formation is not shown for the .S-LC06A/B series, but is 
provided in Table[2] 

In the top panel of Figure [9] we plot the bounce compact- 
ness £2.5 as a function of progenitor ZAMS mass Mzams for a 
range of progenitors from multiple stellar evolutionary stud- 
ies. Even within a given model set, the Mzams -£2.5 mapping 
is highly non-monotonic. At the low end of Mzams covered 
by Figure |9j where mass loss has little influence even in pro- 
genitors of solar metallicity, variations in £2.5 are due predom- 
inantly to particularities in late burning stages, caused, e.g., 
by convective versus radia tive core burning and/ or differences 
in shell burning episodes dWooslev et al.l 120021) . At the high 
ZAMS-mass end, £2.5 is determined by a competition of mass 
loss and rapidity of nuclear-burning evolution. 

The bottom panel of Figure [9] depicts the time to BH for- 
mation ?bh in a failing CCSN as a funct ion of Mzams for 
the sWH07 solar-metallicity progenitors of lWooslev & Hegerl 
(12007b and th e wWHW02 10" 4 solar-metallicity models of 
IWoosley et ail ( 120021) . Models of very low £2.5 that require 
more than 3.5 s to make a BH are omitted. As demonstrated 
in Section |4~4l f BH scales oc (£,2.5)~ 3 ^ 2 and, hence, progen itors 
that form BH the fastest and are (generally, cf. Section 14751 
the hardest to explode are those with high values of £2.5. In 
the low-metallicity mWHW02 series whose progenitors expe- 
rience only minuscule mass loss, BHs form within < 1 s of 
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bounce forMzAMS ^ 30M Q and the high bounce compactness 
£2.5 > 0. 45 m akes a successful shock revival rather unlikely 
(Section 14.5b . Hence, the most likely outcome of core col- 
lapse is BH formation in these progenitors. This may also be 
the case for mWHW02 progenitors in the ZAMS mass range 
from ~ 20 to 25 M Q . The sWH07 progenitors have high £2.5 
and form BHs rapidly only in the Mzams ranges ~ 23-25M 
and ~ 35-45M Q . At higher ZAMS masses, strong O-star 
mass loss leads to an early removal of the hydrogen enve- 
lope. Subsequent mass loss in the W-R phase leads to bare, 
low-mass, low-compactness carbon oxygen cores in the most 
massive progenitors that are unlikely to make BHs. 

4.6.2. Variations with Mass-loss Prescriptions 

Mass loss is key in determining the observational ap- 
pearance of a su ccessful CCSN (e.g., iFilippenkol 119971: 
iSmith etal.ll20lQt) . but, as we have seen in Section 34.6. II 
also has a strong effect on presupernova core structure and, 
thus, on the outcome of core collapse. The det ails of mass 
loss in massive stars are still rather uncertain (e.g. JSmith et al.l 
120101) . and, unfortunately, there are few stellar evolution stud- 
ies that have studied the effects of var iations in mass-loss 
prescriptions. iLimongi & ChieffH ([2006J3 performed such a 
study, adopting two different mass-loss rates for the W-R 
stage of solar-metallicity stars with M > 40 M Q . The 5 LCO6B 
model s are evolved with the W-R mass-loss rates of iLangerl 
(119891) that are similar to those used in the sWH07 set of 
iWooslev & Hegerl ((2007). As depicted in the top panel of 
Figure [5] high-mass 5LCO6B and sWH07 models have sim- 
ilar low £2.5 and most likely do not form BHs but rather ex- 
plode as type-Ibc CCSNe. The models of the sLC06A set 
were ev olved with the lower (fa ctor of ~ 2) W-R mass-loss 
rates of iNugis & Lamersl (120001) . The 5LCO6A 60, 80, and 
120M Q progenitors have much more mass left at the presu- 
pernova stage (M pre _sN ~ 17-30M Q , Figure|2]i and very high 
bounce compactness of £25 ~ 0.6-0.9. In the likely case of 
CCSN failure, a BH forms within ~ 0.5 s with the LS 1 80 EOS 
and within ~ 1 .5 s for all other EOS. 

The above results highlight the sensitivity of outcome pre- 
dictions on mass-loss physics and a more solid understanding 
of this key ingredient will be necessary to robustly connect 
ZAMS masses to the outcome of core collapse for massive 
stars around and above solar metallicity. 

4.7. The Formation of Rotating Black Holes 

Rotation, if sufficiently rapid, alters the CCSN dynamics 
via centrifugal support. This important effect is captured by 
GRID'S 1.5D rotation treatment, albeit approximately. Ini- 
tially, centrifugal support acts to slow the collapse of the in- 
ner core, delaying core bounce. At bounce, lower peak den- 
sities are reached, the hydrodynamic shock forms at a larger 
radius, and its enclosed mass is larger. Conservation of angu- 
lar momentum spins up the core from precollapse angular ve- 
locities that may be of order rads" 1 to O (lOOOrads -1 ) as the 
core, initially with r ~ O(1000 km), collapses to a PNS with 
r ~ O(30km). During the postbounce evolution, the spinning 
PNS is stabilized at lower densities, is less compact, generally 
colder, and has a softer neut rino spectrum than a non-spinning 
counterpart (Ott et al. 2008 and references therein). 

3 See also lLimonei & Chieffl 1 2009) and lMevnet & MaedeJ <2003l) 
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Figure 10. T/\W\ (left) and central density (p c ) (right) during the post- 
bounce evolution of the u40WHW02 model using the LSI 80 EOS and 14 
different initial specific angular momentum profiles. We vary jis,oo from 
to 3.25 in increments of 0.25. For clarity we highlight with solid lines the 
simulations with integer values of y'i6.oo- Lines at T/\W \ = 0.27 and 0.14 are 
added to denote the dynamical and secular rotational instability thresholds. 



4.7.1. Models with Parameterized Rotation 

We investigate the effect of rotation in failing CCSNe by as- 
signing specific angular momentum profiles to the mWHW02 
model set (see Section 0) via Equation (fTTT) . This rotation 
law approximates what is generally foun d in stellar evolu- 
tio n calculations th at account for rotation dHeger et al.ll2.000h 
see lOtt et alJl2006l for comparison plots). The inner iron core 
(~ IMq) is rotating almost uniformly. Outside of this core, 
the angular velocity drops roughly oc r~ 2 . In Table [5] we sum- 
marize key parameters of our rotating model set. Among them 
is r/|W|, the ratio of rotational kinetic energy to gravitational 
binding energy. It is particularly indicative of the dynamical 
relevance of rotation. 

In the right panel of Figure [lOl we plot the central den- 
sity evolution of model m40WHW02 run with the LSI 80 
EOS for 716,00 ran g m g from to 3.25 in increments of 0.25. 
While we choose the m40WHW02 model here, the results are 
generic and apply to all progenitors. Am @ , of Equation (fTTI) . 
is 936 km for this model and the initial central rotation rate is 
1.14 x y'i6.oorads -1 . The nonrotating model takes ~ 433 ms 
to reach bounce and a further ~ 369 ms before the PNS be- 
comes unstable to collapse to a BH with a gravitational mass 
of 2.24M Q . For the 716,00 = 1> 2, and 3 models, respectively, 
the times to bounce are 1 1 ms, 47 ms, and 125 ms greater than 
in the nonrotating case. Their times to BH formation Jbh 
are 12 ms, 52 ms, and 150 ms longer than in the nonrotating 
case. The maximum gravitational PNS masses M g max are 
0.03 M Q , 0.09 M Q , and 0.28M o greater. We find that the time 
to bounce, time to BH formation, and maximum gravitational 
PNS mass increase above the nonrotating values proportional 
to ~ Oi6,oo) 2 - The increase in ?bh is due almost entirely to the 
increase in M g , max , since the accretion rate is not significantly 
affected by rotation. 

The lower temperatures and densities of rotating PNSs lead 
to systematically lower mean neutrino energies and total ra- 
diated energy from the PNS core (ti me-averaged tota l lumi- 
nosities are summarized in Table |5]). iFrver & Hegerl ( 120001) 
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and lOtt et al.l (120081) . who considered similarly rapidly rotat- 
ing models, also see this effect. There is a clear trend toward 
lower L„ with increasing j\6.oo and for a given model and at 
a given time, with increasing 716,00 > l ess gravitational binding 
energy has been carried away by neutrinos and M g is larger. 
Given essentially unaltered accretion rates, one may expect 
earlier PNS collapse and BH formation. This, however, is not 
the case in models run with the LSI 80, LS220, and HShen 
EOS, since the centrifugally increased M g max systematically 
outweighs the increased gravitational mass due to the lowered 
neutrino emission. For these EOS, the time to BH formation 
is delayed by rotation. For models run with the extremely stiff 
LS357 EOS the situation is different. For them, the centrifu- 
gal support provided by rotation is too weak to significantly 
enhance M g max and, hence, BHs form faster with increasing 
716,00 ■ 

In the left panel of Figure QJJ we plot the 7*/|W| evolution 
for the rotating «40WHW02 model series run with the LSI 80 
EOS. During collapse, gravitational binding energy is trans- 
ferred to rotational energy, increasing the value of T/\W\. 
Similar to how the central density overshoots its new equilib- 
rium, also exhibits a local maximum at bounce. Con- 
tinued accretion and contraction of the PNS increases T/\W\ 
throughout the postbounce evolution for all models. Initially 
very rapidly spinning models experience core bounce under 
the strong influence of centrifugal effects, leading to reduced 
compactness and at bounce. These qualitative fea- 
tures are in good agreement with what was found by pre- 
vious extensive parameter studies of rot ating core collapse 
(lOtt et al.l2006tlDimmelmeier et al.l2008l) . Quantitatively, we 
find and summarize in Table[5]that models with 716,00 1 -5 
yield T/\W\ < 0.14 throughout their entire evolution. Mod- 
els with 716,00 ^ 2.25 have T/\W\ > 0.14 during their en- 
tire postbounce evolution. Models that have 716,00 ^ 2.25 
have T/\W\ < 0.27 at all times. Models with j[ 6tO0 > 2.5 
reach r/|W| > 0.27 before BH formation. When considering 
these numbers, it is important to keep in mind that GRID'S 
1.5D approach to rotation has the te ndency to overe stimate 
7*/|W| in rapidly spinning models. lOtt et al.l (120061) found 
model-dependent differences in T/\W\ of O(10%) between 
1.5D and 2D. In addition, GRID'S neutrino leakage scheme 
also tends to lead to somewhat more compact PNS cores and 
consequently higher T/\W\ than would be expected from full 
neutrino transport calculations 

The systematics of T/\W\ depicted by Figure [10] (left) and 
listed in Table [5] albeit only approximate due to GRID'S 
1.5D treatment of rotation, shed interesting light on the po- 
tential role of nonaxisymmetric rotational instabilities during 
the evolution of failing CCSNe. Of course, due to its 1.5D 
nature, GRID cannot track the development of such multi- 
dimensional dynamics. Analytic theory and to some extent 
3D computational modeling have identified multiple instabil- 
ities that may lead to triaxial deformation of PNSs, redistribu- 
tion of angular momentum, and to the radiation of rotational 
energy and angul ar m omentum via gravitational waves (see 
IStergioulasi 120031 and lOttl 120091 for reviews) A global dy- 
namic al instability sets in for T/\W\ > 0.27 (IChandrasekharl 
119691) . leading to a lowest-order m = 2 "bar" deforma- 
tion. Global secular instability, drive n by viscosity or GW 
back-reaction sets in at T /\W\ > 0.14 (IChandrasekharil 1 9701: 
iFriedman & Schutzl [T978). Finally, dynamical shear instabil- 
ities, arising as a result of differential rotation, may lead to 
partial or global nonaxisymmetric deformation at even lower 



value s of T l\W\ (> 0.05: e g ISaiio et all [200l lOtt et al.1 



120071 IScheidegger etaill2008UCorvino et d.l 12010, and ref- 
erences therein). In nature, and in full 3D simulations, these 
instabilities, through gravitational radiation or redistribution 
of angular momentum, will effectively and robustly prevent 
7*/ [W] from surpassing the corresponding T/\W\ threshold. 
The growth times of dynamical instabilities are short, O(ms). 
Secular instabilities grow on ti mescales set by the driving 
process and are typically O(s) dLai & Shapirol fl995). The 
low-T/|W| shear instabilities in PNSs appea r to grow on in- 
termediate timescales o f 0(10-1 00ms) (e.g.. lOtt et~aTll2007t 
IScheideggereraT1l2008l) . 

In Figure [TT] we plot the value of T/\W\ (left panel) and 
the dimensionless spin of the protoblack hole (PBH), «p BH = 
JpBn/(Mg PBH ) (right panel) at the onset of BH formation 
(when a c = 0.3) for the same values of 716,00 used in Fig- 
ure [TO] Assuming that the entire PNS is promptly swallowed 
once the horizon appears, flp BH corresponds to the BH birth 
spirfl. We again show results for model m40WHW02, but for 
all four EOS. The data are also presented in Table|5]for these 
and other models. r/|W| at BH formation scales oc (716 oo) 2 : 
r/|W| PBH is ~ 0.05, ~ 0.1, - 0.2, and - 0.3 at j 16>o0 of - 1, 
~ 1.5, ~ 2.2, and ~ 2.75, respectively. flp BH scales linearly 
with 716.00. reaching a maximally Kerr value of flp BH ~ 1 at 
7i6,oo ~ 2.75. r/|W| PBH and a* m vary little with EOS. 

A disturbing fact depicted by Figure [TTJ is that our 1 .5D 
simulations predict BH birth spins a* > 1 for 716,00 ^ 2.75. 
In Kerr theory, such BHs cannot exist with a horizon. They 
would instead be nak ed singularitie s, violating the cosmic 
censorship conjecture dPenrosell 19691) . However, when com- 
paring right and left panels of Figure [TTJ one notes that all 
models achieving a* > 1 are predicted to reach T /\W\ in ex- 
cess of 0.27. Hence, in nature and in a 3D simulation, these 
PNS will be dynamically nonaxisymmetrically unstable and 
angular momentum redistribution and gravitational radiation 
will limit their r/|W| robustly below ~ 0.27, correspond- 
ing to a* < 0.9. Rotational instabilities at lower values of 
r/|W| may also be relevant. Dynamical shear instabilities 
have timescales significantly less than the time to BH for- 
mation. Secular rotational instabilities may be relevant if the 
true nuclear EOS allows for a large maximum PNS mass as 
more time is neede d to accrete the necessary material to form 
a BH (see Sectio n l4~3l . Large f B e is also possible if £2.5 is 
small (see Section l4~4i i therefore allowing secular instabilities 
to grow. In all rotating models considered here (see Table [5}, 
PNSs stable against the dynamical rotational instability with 
r/|W| < 0.25-0.27 throughout their postbounce evolution 
have < 0.9. Similarly, PNSs with T/\W\ < 0.14-0.16, 
the threshold for secular instability, have flp BH < 0.6-0.7. 
If low-r/|W| instabilities are effective at limiting T/\W\ in 
PNSs on short timescales, nascent BH spins may be limited 
to low values (a* < 0.4 for a r/|W| instability threshold of 
- 0.05). 

4.7.2. Rotating Progenitors and the Connection to long GRBs 

The rotation law of Equation (fTTT > qualitatively follows the 
predicted angular velocity distribution in the inner ~ 1 - 3M 
of presupernova models evolved with rotation. However, 

4 Note that this may not necessarily be what happens. Outer PNS material 
may become centrif ugally supported, falling into the nascent BH only on an 
accretion timescale (Duez et al. 2006). 
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Note. — The u series of presupernova models in this table are taken from the Zq = 10~ 4 model set of Woosley et al. 1 2002). We imposed a rotat i on law via 
equation [T7I the value of y'u x is giv en in the model na me following the letter J. The m350C presupernova model is taken from Wooslev & Heser ( 2 006T) and 
both the E20 and E25 models are from Heser et al. 1 2000), these model are evolved with rotation. In simulations where a BH did not form within 3.5 s we give, in 
parenthesis, the values at this time. 
a Initial central angular velocity of star. 
b Total angular momentum of star. 
c Initial T/\W\ of the star. 
d Tj \W\ of the star at bounce. 

e Angular momentum of protoblack hole (PBH) when a c = 0.3. 
f T/\W\ of the star when a c = 0.3. 

g Dimensionless spin of the PBH when a c = 0.3. Unphysical values for a BH are shown in braces [■ ■ ■]. 
h Total neutrino luminosity averaged over postbounce time. 



17 



0.5 



0.4 



_0.3 
^0.2 



0.1 



1 1 1 1 1 1 1 1 1 1 1 1 1 

«40WHW02 g . 

• LS180 

■ LS220 *- 

~ LJJ / J 

* HShen * ■ 

■ 

"T/IWI = 0.27 - * 


f 

i 

§...- 

1 

■ 

■ 


A 

! 

T/IWI = 0.14* 


• 

■ 

m 

■ 


1 

1 

■ 

■ 

.-•,1,1,1,1,1,, 


■ 

• 

" ■ 

,,1,1,1,1,1,1, 



1.25 



0.75 



0.5 



-0.25 



ni-»j 4 o 

0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 3.5 

j 16oo [10 16 cm 2 s 1 ] 

Figure 11. Left: 7"/|W| for a range of initial 716, 00 and EOS for the 
u40WHW02 progenitor. We denote the value of T/\W\ thresholds for the 
dynamical rotational instability, r/|W|;j yn = 0.27, and the secular instability, 
r/|W| S cc = 0.14. Right: dimensionless spin parameter flp BH for the PNS at 
the last stable configuration prior to collapse to a BH. Note that £'p BH > 1 
is generally allowed by GR but a BH must have a* < 1. PNSs that could 
reach cip BH > 1 are nonaxisymmetrically unstable and will be limited to cip BH 
below 1 . For the uWHW02 model, the initial central rotation rate is given as 
Sl ( . = 1.141 X ;'i6, ooi'ads -1 . 



Equation (fTTT i asymptotes to constant specific angular mo- 
mentum j and cannot capture jumps and secular increase of j 
in overlying mass shells (e.g. jHeger etaill2000l) . 

Here we consider three different supernova progenitors 
evolved with rotation that have the potential of forming BHs 
soon after bounce. Models E20 and E25 are rapidly rotat- 
ing unmagnetized solar-metallicit y presupernova mo dels of a 
20 and 25 Mq ZAMS stars from iHeger et all (120001) . Model 
m350C o flWooslev & He gerl (120061) withM 7A Ms = 35 M Q has 
10 % solar-metallicity, reduced mass loss, and magnetic fields. 
These presupernova models have initial central angular ve- 
locities of ~ 3.13, ~ 1.83, and ~ 1.98rads _1 and values of 
£2.5 of ~ 0.319, ~ 0.294, and ~ 0.456 for the E20, E25, and 
m350C models, respectively. Due to their moderate £2.5, we 
perform collapse simulations of models E20 and E25 with the 
LSI 80 and LS220 EOS. Model m350C is calculated with all 
four EOS. The progenitor characteristics are summarized in 
Tables Q]and[5] 

For a given EOS, due to very similar £2.5.. the evolutions 
of E20 and E25 are alike. They form BHs in ~ 1.1 s with 
the LSI 80 EOS and in - 1.7s with the LS220 EOS. Model 
E20 is more rapidly spinning. Its r/|W| peaks at ^0.18 at 
bounce and settles down to a nearly constant value of ~ 0.15 
throughout the postbounce evolution. It's flp BH is ~ 0.7 for 
both the LS 1 80 and LS220 EOS . Model E25 reaches T/\W\ ~ 



0.06 at bounce and < 
for both EOS. 



■ 0.065 at BH formation with a 



PBH 



-0.45 



The core of the m350C model is sufficiently compact to 
form a BH soon after bounce if no explosion is launche d (e.g., 
via magneto-rotational explosion; iDessart et al . 1 120081) . The 
nascent BH forms at a time of ~ 0.75 s, ~ 0.97 s, ~ 2.19s, 
and - 1.91 s for the LSI 80, LS220, LS375, and HShen EOS, 
respectively. The initial gravitational (baryonic) BH masses 
are - 2.22 (~ 2.33)M Q , - 2.29 (~ 2.43)M , ~ 2.69 (~ 
3.O1)M , and ~ 2.60 (~ 2.84)M Q for the LS180, LS220, 



LS375, and HShen EOS, respectively. The BHs are modestly- 
rapidly spinning with flp BH of <~ 0.61, ~ 0.58, ~ 0.54, and 
~ 0.55 for the LSI 80, LS220, LS375 and HShen EOS, re- 
spectively. For all EOS, the PNS, during the accretion phase, 
has a modest T/ \ W\ of < 0.12. 

Once a BH is formed, material from the stellar mantle 
will continue to accrete at high rates. Accretion will only 
be slowed once material with sufficiently high specific an- 
gular momentum reaches small radii and becomes centrifu- 
gally supported, forming an accretion disk. This is the crucial 
prerequisite for t he collapsar sce nario for long GRB central 
engines to work (IWoosleylll993l) . Models E20 and E25 lost 
much of their initial mass and angular momentum during their 
evolution to the presupernova stage and, therefore, there is too 
little angular momentum in their outer regions to allow for a 
long-term accretion disk. 

The situation is different for the m350C progenitor. Its par- 
ticular evolution prevented dramatic loss of mass and angular 
momentum while keeping its envelope radius small. In Fig- 
ure [12J we show the specific angular momentum distribution 
of the m350C progenitor as a function of enclosed baryonic 
mass. We also include graphs of the yisco, the specific angular 
momentum required for a stable orbit at the innermost stable 
cir cular orbit (ISCO) o f a Kerr hole with mass M and spin 
a* dBardeen et alj|1972h . Curves for a* = (Schwarzschild), 
a* = 1 (maximally Kerr), and a* = J(M)/M 2 are shown. We 
also plot the value of a* that a BH of baryonic mass M formed 
from the m350C progenitor would have. Figure [12] is inde- 
pendent of the detailed collapse evolution, assuming that no 
angular momentum is radiated by neutrinos and/or gravita- 
tional waves, or ejected. However, we note that due to emis- 
sion of neutrinos before BH formation, the enclosed gravita- 
tional mass (entering into the calculation of a*) will be smaller 
by up to ~ 0.2-0.4Mq than the baryonic mass given in the 
figure. This leads to slightly underpredicted values of a* for 
small M. Since relatively little energy is emitted in neutrinos 
after BH formation, the relative discrepancy between gravita- 
tional and baryonic mass decreases with growing BH mass. 

Figure[T2]can be interpreted as follows. If the CCSN mech- 
anism fails to reenergize the shock and the PNS collapses to a 
BH of mass M, then its initial angular momentum J and spin 
a* will be set by the enclosed angular momentum and gravita- 
tional mass. Initially, hyperaccretion will increase both J and 
M, but a* may increase or decrease, depending on the angular 
momentum of the accreted matter. Accretion will slow down 
and a disk will form once infalling material has specific an- 
gular momentum j greater than j'isco- In model m350C, this 
occurs between a BH mass coordinate of ~ 7.6M (for a mass 
element with equatorial j) and ~ 10. 1 M (for a mass element 
with angle-averaged j). These points are marked in Figure [T2l 
with a (-^r) and (■), respectively. Using accretion simula- 
tions with GRID setup to include an inflow inner boundary 
condition we find that with the m350C model, the accretion 
time for 7.6M and 10.1M S is ~ 10.2s and — 14.3s from 
the onset of collapse, respectively. These times are roughly 
twice the free fall time, since material in out er regions in hy - 
drostatic equilibrium for a sound travel time dBurrowsll 19861) . 
Once the disk has formed, accretion will continue via pro- 
cesses that will transport angular momentum out an d mass in. 
A col l apsar central engine may begin its operation (IWoosleyl 
119931: iMacFadven & Wooslevl [1999b with a central BH of 
M - 8M , a* ~ 0.75, and an ISCO radius of - 40km. 
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Baryonic Mass [M Q ] 

Figure 12. Specific angular momentum (solid, left ordinate) and dimen- 
sionless spin (dot-dot-dashed, right ord inate) for the GRB progenitor model 
m350C from Woosley & Hegei 12006) as a function of enclosed baryonic 
mass. The thin solid line is the angular momentum at the equator; the 
thick solid line is the angle-averaged angular momentum. Also shown, dash- 
dotted, dash-dash-dotted, and dashed, is the specific angular momentum for a 
mass element at the innermost stable orbit assuming a background spacetime 
being Schwarzschild, maximally Kerr, and Kerr with a* =J(M)/M 2 , respec- 
tively. The horizontal line denotes a* = 1, (+) and (■) denote the mass 
coordinate where the equatorial and angle-average specific angular momenta 
exceed Jisco respectively. 



5. SUMMARY AND CONCLUSIONS 

We have performed an extensive study of BH formation in 
failing CCSNe with the open-source 1.5D GR code GRID, 
making the simplifying assumptions of spherical symmetry 
and of a neutrino leakage scheme rather than full Boltzmann 
transport. We have performed more than 700 collapse simu- 
lations with over 100 unique progenitor models, probing sys- 
tematically the many-dimensional parameter space that deter- 
mines the outcome of stellar collapse in single massive stars. 
Specifically, we have studied and established the systematic 
dependence of CCSN failure and BH formation on progeni- 
tor compactness, precollapse rotational setup, neutrino heat- 
ing efficiency, and nuclear EOS. 

To first approximation, the evolution of any core collapse 
event proceeds as follows. Core collapse ensues in a given 
presupernova star, collapse is halted when the inner core of 
~ 0.5-0.7M© reaches nuclear density. A shock is formed, 
propagates outward initially in M and r, but soon stalls. As- 
suming the CCSN mechanism, whatever its precise nature 
may be, fails, we can robustly predict the time it takes to BH 
formation for a given nuclear EOS (scaling with EOS stiff- 
ness) based on a single parameter, the progenitor bounce com- 

pactness £2.5 (?bh <x £2.5 )■ Using the same parameter, for 
a given EOS, we can predict the maximum mass of the PNS 
at collapse and its thermal enhancement (10%-25%) over the 
cold NS mass, due, as we have shown for the first time by 
detailed comparison with exact TOV solutions, primarily to 
thermal pressure support in the outer PNS core. 

In an attempt to more quantitatively understand which stars 
explode and which do not, assuming the neutrino mecha- 
nism is responsible for the majority of CCSN explosions, we 
have turned the knobs on GRID'S neutrino heating scheme, 
experimentally, to first order, establishing the neutrino heat- 



ing efficiency needed to explode a progenitor with given £2.5- 
Neglecting the potentially highly relevant effects of multi- 
dimensional dynamics and assuming an EOS of intermediate 
stiffness (the LS220 EOS), we predict that progenitors with 
bounce compactness £2.5 si 0.45 most likely form BHs with- 
out explosion. This prediction, in itself, without connection to 
ZAMS conditions through stellar evolution, is of limited util- 
ity. Using the whole set of progenitor data made available to 
us by stellar evolution groups, we attempt the former in Fig- 
ure Qj] We plot the mapping between ZAMS mass and out- 
come of core collapse, reduced to explosion or no explosion 
and BH formation, neglecting completely the possibility of 
BH formation due to fallback/cooling/phase transitions after 
a launched explosion. The case is clear cut at low metallic- 
ity where mass loss has negligible effect on the mapping be- 
tween ZAMS conditions and core collapse outcome. Using a 
Salpeter initial mass function (IMF; a = 2.35, M min = 8.0M©, 
and M max = 150. 0M Q ) we estimate that ~15% of all progen- 
itors form BHs without explosion. At (around) solar metal- 
licity, the precise way of prescribing mass loss in stellar evo- 
lution has tremendous consequences on the mapping between 
ZAMS mass and core collapse outcome. Depending on the 
particular mass-loss prescription, we predict a BH fraction 
of 0%-7% for solar-metallicity stars. This makes mass loss 
the single most important unknown parameter in connecting 
ZA MS conditions to core collapse outcome (in agreement 
with lSmith e"taill2010l) . 

Rapid rotation, which may be present in a significant subset 
of massive stars, generally increases the maximum PNS mass 
by centrifugal support and delays BH formation. Assuming 
(quite likely) uniform rotation of the PNS core, the increase in 
maximum PNS mass due to centrifugal support in the range 
of rotation rates explored is ~ 5%-10%. In the basic neu- 
trino mechanism, rotation leads to a lower sum of v e and v e 
luminosities and lower mean energies for all neutrino types. 
This is detrimental for explosion in 1 , 5D (and perhaps even 
in 2.5D) desp ite centrifugal support dFrver & Hegerl 120001 : 
lOtt et al.ll2008l) . A larger fraction of massive stars may form 
BHs with (moderate) rotation than without. Left out of this 
picture are potential magnetohydrodynamics c ontributions to 
the exp losion mechanism and energetics (cf. iBurrows et al.l 
I2007bl> . 

Of particular interest to both formal relativity theory and 
astrophysics is the range of potential birth spins of BHs. 
Our results quite strikingly suggest that the rotation rate of 
the maximum-mass PNS and, hence, the spin of the nascent 
BH, will be limited to values of a* below < 0.9 by likely 
nonaxisymmetric dynamics. If true and confirmed by multi- 
dimensional simulations, 3D rotational instabilities may be a 
cosmic censor preventing naked singularities from forming in 
stellar collapse. 

Rotation and the associated angular momentu m are key 
ingre dients in the collapsar scenario for GRBs dWooslevI 
1993). As part of this study, we have performed the first 
BH formation study wi th the m350C GRB progenitor of 
IWooslev & Hegerl (120061) . Using the LS220 EOS, we predict 
an initial BH mass of ~ 2.29 M Q and a* of ~ 0.58. Assuming 
that the GRB engine cannot operate until a Keplerian disk has 
formed, there will be a delay of ~ 10 s between BH forma- 
tion and GRB engine ignition at a BH mass of ~ 8M Q and 
a* -0.75. 

Finally, we re-emphasize that the goal of this study was 
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Figure 13. Outcome of core collapse as a function of ZAMS mass of single nonrotating massive stars, assu ming that for moderately stiff nuclear EOS (e.g., 
LS 1 80/LS220), neutrino-driven explosions can be launched up to a bounce compactness §2.5 5; 0.45 (cf. Section RTBl . Other potential explosion mechanisms are 
neglected. We consid er only explosion and BH formation w ithout explosion as outcomes and neglect other scenarios, including post-explosion BH formation 
via fallback accretion iZhang et al. 2008; Dessart et al. 2010), cooling or nuclear phase transitions. Shown are results for a range of model sets and metallicities 
(see Section[3). Very low met allicity stars w ith ZAMS masses above ~ 30Mq robustly form a BH without explosion. At higher metallicity, uncertainties in the 
physics of mass loss (e.g., Smith et al. 2010) make robust predictions difficult. This is reflected in the rather dramatic disagreement of the four solar-metallicity 
progenitor model sets that we include. The "BH fractions" stated at the right edge of the plot denote the fraction of massive stars with M > &Mq that form BHs. 
They are obtained by convolution with a Salpeter IMF under the assumption that stars with 8Mq < M < 14Mq explode robustly. 



not to yield accurate predictions about the outcome of core 
collapse in any individual progenitor. Rather, we have stud- 
ied and established overall trends with progenitor parameters. 
We have made simplifications and approximations, and have 
omitted a broad range of potentially relevant physics. The 
most important of the latter may well be multi-dimensional 
dynamics and their effect on the CCSN explosion mechanism 
and on the associated failure rate of CCSNe. 

Future work may be directed toward studying the system- 
atics of BH formation in the post-explosion phase via fall- 
back accretion, PNS cooling, or EOS phase transitions. Our 
current neutrino treatment must be upgraded for more quan- 
titatively accurate simulations and neutrino signature predic- 
tions. Ultimately, multi-dimensional GR simulations of suc- 
cessful and failing CCSNe will be necessary to study the 
multi-dimensional dynamics left out here and for making truly 
robust predictions of the outcome of stellar collapse for any 
given set of initial conditions. 
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